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Abstract 

Methods based on Wilson's renormalization group have been successfully applied in the context of nuclear physics 
to analyze the scale dependence of effective nucleon-nucleon (A^A^) potentials, as well as to consistently integrate out 
the high-momentum components of phenomenological high-precision A^A^ potentials in order to derive phase-shift 
equivalent softer forms, the so called Vunv-k potentials. An alternative renormalization group approach that has been 
applied in this context is the Similarity Renormalization Group (SRG), which is based on a series of continuous 
unitary transformations that evolve hamiltonians with a cutoff on energy differences. In this work we study the SRG 
evolution of a leading order (LO) chiral effective A^A^ potential in the '5o channel derived within the framework of the 
Subtracted Kernel Method (SKM), a renormalization scheme based on a subtracted scattering equation. 

Keywords: similarity renormaUzation group, chiral effective field theory, nucleon-nucleon interaction. 



1. Introduction 

Chiral Effective Field Theory (ChEFT) offers a promising framework for a systematic and model-independent 
description of hadronic and nuclear systems and their interactions in the low-energy regime that is consistent with 
Quantum Chromodynamics (QCD), the fundamental theory of strong interactions (for detailed reviews, see e.g. Refs. 
||T}{4|). The construction of such a framework began with the advent of Chiral Perturbation Theory (ChPT), proposed 
by Weinberg |5| and successfully implemented in calculations of various low-energy processes forpion-pion (nn) and 
pion-nucleon {nN) systems ll6 Wn . As shown by Weinberg in his seminal paper |5|, starting from the most general 
lagrangian for pions and nucleons consistent with the underlying symmetries of QCD, in particular the approximate 
and spontaneously broken S U(2\ x S U(2)r chiral symmetry, a systematic low-energy expansion can be derived 
for the 5 -matrix in terms of powers v of the ratio (Q/A^), where Q stands for a low-momentum scale associated 
with external pion and nucleon momenta or the pion mass and A^ ~ 1 GeV is the chiral symmetry breaking scale. 
For purely-pionic and single-nucleon systems the interactions become weak in the low-energy regime, due to the 
Goldstone boson nature of pions, ensuring that the chiral expansion of the 5 -matrix converges. Such a perturbative 
expansion corresponds to an infinite number of Feynman diagrams which can be organized according to a power 
counting scheme based on naive dimensional analysis. At any given order v in the chiral expansion, there is only 
a finite number of diagrams. This power counting scheme, known as Weinberg's power counting, is an essential 
ingredient for systematic and controlled calculations of low-energy observables in the ChPT approach, which allows 
to select the relevant diagrams contributing to the 5 -matrix up to a given order and to estimate the truncation errors 
involved in the calculations. 

The success of the ChPT approach in the nn and nN sectors led Weinberg lfT2ffl4l to propose its extension to 
describe the nuclear forces. However, there is a fundamental difficulty with the direct application of ChPT methods in 
the few-nucleon sector, where the interactions are strong in the low-energy regime. The existence of shallow nuclear 
bound states (e.g. the deuteron) and large scattering lengths in 5 -wave channels clearly indicates the non-perturbative 
nature of the nucleon-nucleon (A^A^) interaction at low-energies. Thus, a purely perturbative treatment of few-nucleon 
systems is expected to fail. 
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Weinberg pointed out that the non-perturbative dynamics of few-nucleon systems is due to the strong enhancement 
of the S -matrix arising from purely nucleonic intermediate states (infrared enhancement) and then outHned a two-step 
strategy to perform non-perturbative EFT calculations [12. , 131. In the first step, the power counting scheme of ChPT 
(Weinberg's power counting) is applied to an effective nuclear potential rather than directly to the full S -matrix. For 
a «-nucleon process, the effective potential is defined as the sum of all possible n-nucleon time-ordered perturbation 
theory diagrams without purely nucleonic intermediate states, the so called irreducible diagrams, obtained from the 
chiral expansion of the lagrangian for jttt and ttN interactions supplemented by a lagrangian for n-nucleon contact 
interactions constrained only by isospin symmetry. In the second step, the effective potential truncated at a given 
order in the chiral expansion is inserted into the Lippmann-Schwinger (LS) or the Schrodinger equation and then 
iterated to all orders to obtain the full 5 -matrix, thus generating the non-perturbative effects. 

The effective nuclear potential includes long-range contributions from pion exchange interactions and short-range 
contributions parametrized by nucleon contact interactions with an increasing number of derivatives. Weinberg's 
power counting, used to organize the chiral expansion, naturally implies in the observed hierarchy of the nuclear 
forces {14]. At leading order (LO, v - 0), the effective nuclear potential consists of the well known A^A^ one-pion 
exchange potential (OPEP) plus two non-derivative A^A^ contact interactions. Contributions from multi-pion exchange 
interactions start at next-to-leading order (NLO, v - 2) and continue through all higher orders. Derivative contact 
interactions also start to contribute at NLO and due to parity appear only in even orders. According to Weinberg's 
power counting, contributions from three-nucleon (3A^) and four-nucleon (4A^) forces should start at NLO. However, 
due to additional suppressions and cancelations (both in energy-dependent and energy-independent formulations) the 
first non-vanishing contributions from 3A^ and 4N interactions occur, respectively, at next-to-next-to-leading order 
(NNLO, V = 3) and next-to-next-to-next-to-leading order (N^LO, v = 4) 12, ,4 1. Both pion exchange and contact 
interaction terms usually grow with momenta and thus can lead to ultraviolet divergences when the effective potential 
is iterated to all orders in the LS or the Schrodinger equation, requiring the use of a non-perturbative regularization 
and renormalization procedure in order to obtain well-defined finite solutions. Furthermore, multi-pion exchange 
interactions involve ultraviolet divergent loop integrals which must be consistently regularized and renormalized. 

The standard procedure for the non-perturbative renormalization of the A^A^ interaction in the context of Weinberg's 
approach to ChEFT consists of two steps |l2l|4l[T5J. The first step is to solve the regularized LS equation for the 
scattering amplitude with the A^A^ potential truncated at a given order in the chiral expansion. The most common 
scheme used to regularize the LS equation is to introduce a sharp or smooth regularizing function that suppresses the 
contributions from the potential matrix elements for momenta larger than a given cutoff scale, which separates high- 
energy/short-distance scales and low-energy/long-distance scales, thus eliminating the ultraviolet divergences in the 
momentum integrals. The second step is to determine the strengths of the contact interactions, the so called low-energy 
constants (EEC's), by fitting a set of low-energy scattering data. Once the LECs are fixed for a given cutoff, the LS 
equation can be solved to evaluate other observables. Such a procedure, motivated by Wilson's renormalization group 
lfT6l - [T8l . relies on the fundamental premise of EFT's (and essentially all renormalization techniques) that physics 
at low-energy/long-distance scales is insensitive with respect to the details of the dynamics at high-energy/short- 
distance scales ifTSll . i.e. the relevant high-energy/short-distance effects for describing the low-energy observables can 
be captured in the cutoff-dependent EEC's. The A^A^ interaction can be considered properly renormalized when the 
calculated observables are independent of the cutoff scale within the range of validity of the ChEFT or involves a small 
residual cutoff dependence due to the truncation of the chiral expansion. In the language of Wilson's renormalization 
group, this means that the EEC's must run with the cutoff scale in such a way that the scattering amplitude becomes 
(approximately) renormalization group invariant. 

Since its pioneering implementation, by Ordonez and van Kolck ||T9l , Weinberg's program has been extensively 
explored to analyze the nuclear forces. Effective A^A^ and 3A^ potentials have been derived to several orders in the 
chiral expansion and successfully applied in many calculations of nuclear systems [20 -53 1 . At N^LO, effective 
nuclear potentials have been constructed Il39ll42l which provide a remarkably accurate description of low -energy A^A^ 
scattering phase-shift data and deuteron properties, comparable to that obtained with phenomenological high-precision 
potentials like the Argonne V18 [54J and the Nijmegen [55 1. 

In spite of its phenomenological success, conceptual problems have been raised about Weinberg's approach which 
are related to the formal inconsistency between the naive power counting scheme (Weinberg's power counting) and the 
proper renormalization of the A^A^ interaction in the non-perturbative regime. Such inconsistency arises from the non- 
renormalizability of ChEFT in the usual sense: the iteration of the A^A^ potential truncated at a given order in the chiral 
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expansion (using the LS or the Schrodinger equation) generates ultraviolet divergences, and hence regularization scale 
dependencies, which cannot be absorbed by renormalizing the LEC's of the contact interactions included at that same 
order. Thus, upon cutoff regularization, the limit of infinite cutoff scale cannot be taken while keeping the results 
cutoff-independent. It is important to observe that in successful calculations of A^A^ systems following Weinberg's 
approach, such as those described in Refs. [^71 1221 [33H42ll . cutoff scales ~ 400 - 800 GeV are typically chosen and 
nearly cutoff-independent stable results are obtained provided the cutoff is varied only over a narrow range of values, 
in agreement with the prescription advocated by Lepage ifTSll that the cutoff scale should be taken below A^, ~ 1 GeV. 

The non-perturbative renormalization of the A^A^ interaction in ChEFT has been intensively investigated by many 
authors |,56-, - ,l 14j . with alternative approaches and power counting schemes being proposed, generating a great deal of 
discussion and controversy. In particular, methods based on Wilson's renormalization group have been successfully 
applied to analyze the scale dependence of effective A^A^ potentials, both in momentum ||90H99l and in coordinate space 
mOllfTOSJ . providing a better understanding of the interplay between the power counting and the renormalization of 
the NN interaction in the non-perturbative regime. Although much progress has been made, the construction of a 
consistent scheme for the non-perturbative renormalization of the A^A^ interaction still remains an open problem [ 11 4J . 

Renormalization group methods have also been used to consistently integrate out the high-momentum components 
of both ChEFT and high-precision A^A^ potentials in order to derive phase-shift equivalent softer forms, the so called 
Viow-k potentials 1 1 15 -1241. An alternative renormalization group approach that has been applied in this context is the 
Similarity Renormalization Group (SRG), developed by Glazek and Wilson M125II126J and independently by Wegner 
11271 . The SRG is a renormalization method based on a series of continuous unitary transformations that evolve 
hamiltonians with a cutoff on energy differences. Such transformations are the group elements that give the method 
its name. Viewing the hamiltonian as a matrix in a given basis, the similarity transformations suppress off-diagonal 
matrix elements as the cutoff is lowered, forcing the hamiltonian towards a band-diagonal form. An important feature 
of the SRG is that all other operators are consistently evolved by the same unitary transformation. Recently, the 
SRG has been applied to evolve several A^A^ potentials to equivalent softer forms, effectively decoupling low-energy 
observables from high-energy degrees of freedom II128H1321 . As shown in these references, such a decoupling is 
universal and leads to more perturbative A^A^ potentials, improving variational convergence and greatly simplifying 
calculations in nuclear few- and many-body problems. A recent review on the application of renormalization group 
methods to nuclear forces and nuclear structure calculations can be found in Ref. II133I . 

In this work we study the SRG evolution of chiral effective A^A^ potentials derived within the framework of the 
Subtracted Kernel Method (SKM) 111341 4139 1. The SKM is a renormalization scheme designed to treat interactions 
containing regular and/or singular contact terms which is based on a subtracted scattering equation. Instead of using a 
cutoff regularizing function, in the SKM the scattering equation is regularized by performing subtractions in the kernel 
at a given energy scale, while keeping the original interaction intact. An advantage of the SKM approach is that it 
can be recursively extended to any derivative order of the contact interactions through an iterative process involving 
multiple subtractions. A similar approach based on subtractive renormalization of the LS equation is described in 
Refs. 1 140 -142.1 . although there a sharp momentum cutoff is also introduced to regularize the momentum integrals. 

This paper is organized as follows. In Section |2] we review the basics of the SRG formaUsm. In Section |3] we 
discuss a simple example in order to illustrate the application of the SRG approach: the Schrodinger equation with a 
two-dimensional Dirac-delta contact interaction. In Section[4]we review the main aspects involved in the application 
of the SRG approach to the A^A^ interaction and present the results obtained for the SRG evolution of the Nijmegen 
potential in the '5o channel. From the SRG evolved potentials we calculate the corresponding phase-shifts and use 
an exponential regularizing function II 1291 11301 to verify the explicit decoupling of low-energy observables from 
the dynamics of the high-energy degrees of freedom. In Section |5] we describe the SKM formalism and discuss 
its recursive extension to any derivative order of the contact interactions. We illustrate the application of the SKM 
formalism by presenting two simple examples, namely the A^A^ interaction in LO "pionless" EFT, which consists 
of just a pure Dirac-delta contact interaction, and in NLO "pionless" EFT, which consists of a Dirac-delta contact 
interaction plus a second order derivative contact interaction. In Section|6]we apply the SRG transformation to evolve 
an effective A^A^ potential obtained by implementing the SKM approach to renormalize the A^A^ interaction. First, we 
present a detailed and systematic analysis of the SKM approach applied to the A^A^ interaction in LO ChEFT. Then, we 
evolve the SKM-LO ChEFT potential in the '5o channel through the SRG transformation, calculate the corresponding 
phase-shifts and analyze the decoupling pattern between low- and high-momentum components. Finally, in Section [7] 
we summarize the results and present our main conclusions. 
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2. Similarity Renormalization Group Formalism 

2.1. Glazek-Wilson Formulation 

The general formulation of the SRG was developed by Glazek and Wilson 01251 11261 in the context of light- 
front hamiltonian field theory (HLFFT), aiming to obtain effective hamiltonians in which the couplings between high- 
and low-energy states are eliminated, while avoiding the appearance of artificial divergences in the form of small 
energy denominators in perturbative expansions. The method is implemented by continuous unitary transformations 
designed to replace the effects of the removed couplings by scale dependent effective interactions with no small energy 
denominators, ensuring that the effective hamiltonians produce no ultraviolet divergences. It is important to note that 
the transformations do not remove any degrees of freedom, but only integrate out couplings between states with energy 
differences larger than a given scale. 

Consider a system described by a canonical hamiltonian written in the form H - h + V , where h is the free 
hamiltonian and V is an interaction. In general, the hamiltonian contains direct couplings between states from all 
energy scales which can be the source of ultraviolet divergences. In the following discussion, we use the basis of 
eigenstates of the free hamiltonian, h\ i >= e,| / >. 

Starting from the canonical hamiltonian H we define a bare hamiltonian H,i^ regularized by a large initial cutoff 
Af) on energy differences at the interaction vertices. Such a cutoff is called similarity cutoff and here has dimensions 
of energy. The bare hamiltonian can then be written in the form 

= h + VA, = h+f,„{y + Hl), (1) 

where fx„ is called similarity function and H"^ are counterterms to be determined through the renormalization process 
in order to remove dependence in observables upon the cutoff Aq. The similarity function /i„ is defined to regularize 
the hamiltonian by suppressing matrix elements between states with large energy differences. Usually, the similarity 
function is chosen to be a smooth function of the cutoff, such that 

(04 ^ 1, for |e; - ej\ « Aq , 

(ii)f,^ ^ 0, for |e,. - ej\ » Ao . (2) 

A simpler choice for the similarity function is a step function d(Ao - - e,|). Although useful in analytic calculations, 
such a choice can lead to pathologies [147]. 

The next step is to define a unitary transformation U{A, Aq) that acts on the bare hamiltonian Hi^ and changes the 
similarity cutoff down to a scale A, producing a renormalized hamiltonian, 

H, = U(A,Ao)H,,UHa,Ao) ■ (3) 

The unitarity condition satisfied by U(A, Aq) is given by 

UiA,Ao)uHA,Ao) = U\A,A^,)U{A,Ao)=l . (4) 

The renormalized hamiltonian H\ is driven towards a band-diagonal form as the similarity cutoff A is lowered and 
can be written in the form 

H, = h+f,VA, (5) 

where is the similarity function at the cutoff scale A and Va is defined as the reduced interaction (i.e., the effective 
interaction at the cutoff scale A with the similarity function factored out). 

The similarity transformation U (A, Aq) can be defined in terms of an anti-hermitian operator 77^ which generates 
infinitesimal changes of the cutoff scale A, 

U(A,Aq) = Texp{^j^\,,dA'^ , (6) 

where T puts the operators in order of increasing cutoff scale A' from left to right. 
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Taking the derivative of Eq. (j6| and using both the unitarity condition given in Eq. Q and its derivative we get 
the relation 

tj^^U(A,Aq) — = — U'(A,Ao) ^ -T]\ . (7) 

dA dA 

Thus, taking the derivative of Eq. ([3]) and using the unitarity condition combined with Eq. (j?]) we obtain a first-order 
differential equation for the evolution of the hamiltonian, 

^ = (8) 

which is to be solved with the boundary condition H^l^^^^ = 

The generator tj^ is defined by specifying constraints on the change of the operators h and Va with the similarity 
cutoff A. One possible choice is to demand that h is independent of A and that the renormalized hamiltonian does not 
contain any small energy denominators. Such constraints are defined by the conditions 

dh dV 1 

dA-='^ ^MV.,,]. (9) 

Since the similarity transformation is unitary, both and //^ produce the same spectra for the observables. 
Moreover, when the transformation is implemented in an exact form (i.e., not using a perturbative expansion) the 
observables calculated using the renormalized hamiltonian //^ are independent of the similarity cutoff A. As pointed 
before, the counterterms H"^ are adjusted such that the observables also become independent of A^ in the limit Aq ^ oo. 

An important feature of the SRG is that all operators are consistently renormalized by using the same unitary 
transformation, i.e. any given operator O evolves with the similarity cutoff scale A according to 

'^^[o.M- (10) 

dA 

2.2. Wegner Formulation 

In the applications described in this work, we employ the formulation for the SRG developed by Wegner F127I 
and applied within the context of many-body problems in condensed matter physics 1 14811 . Wegner's formulation is 
based on a non-perturbative flow equation that governs the unitary evolution of a hamiltonian with a flow parameter s 
that ranges from to oo, 

^^[lls,Hs], (11) 

ds 

which is to be solved with the boundary condition //s|„,^ = H,,,. 

The flow equation Eq. ( [TT] i is analogous to Eq.(|8]l, but the specific form rjs = [Gs, Hs] is chosen for the anti- 
hermitian operator that generates the unitary transformation, which gives 

^ = [[G„ . (12) 
ds 

The operator G, defines the generator 77, and so specifies the flow of the hamiltonian. Wegner's choice in the 
original formulation is the full diagonal part of the hamiltonian in a given basis, G , - diag{Hs). A simpler choice is to 
use the free hamiltonian, Gs = h. Such choices for the generator t], result in a flow parameter s which has dimensions 
of (energy)"^ and correspond in the Glazek- Wilson formulation to the choice of a gaussian similarity function with 
uniform width A. In terms of the similarity cutoff A, the flow parameter is then given by the relation s - A^^. Although 
both these choices for Gj are able to force the hamiltonian towards a band-diagonal form as the flow parameter s 
increases (or as the cutoff A decreases), it has been shown by Glazek and Perry 1 149 1 that using Gs - h can produce 
divergences in theories with bound-states. Such an effect is related to the renormalization group limit cycle behavior 
B150II1511 and happens when A = 1/ yfs approaches a scale at which a bound-state emerges, thus limiting how far the 
SRG transformation can be implemented. On the other hand, by using Gs - diag{Hs) the limit cycle behavior can be 
properly handled and the transformation is guaranteed to converge. 
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This is an important issue to be considered in the applications of the SRG to nuclear physics. As has been recently 
shown in several papers M1281413 11. the SRG transformation using G , = h can be applied to evolve the A^A^ interaction 
with no convergence problems as long as A is kept larger than the bound-state scales. On the other hand, in cases with 
limit cycle behavior such as the three-nucleon problem II152H1551 . the SRG transformation must be implemented by 
including interactions in G , |T49|. 

The consistent renormalization group evolution of many-body nuclear forces is a crucial and challenging issue 
for nuclear structure calculations using low-momentum interactions II133I . As shown in recent works, the SRG flow 
equations provide an approach to evolve three-body nuclear forces directly, i.e. without having to solve the full three- 
body problem as required in the V/ok-jI: approach. Preliminary studies using simple models 11561 1157|| have been 
carried out aiming to set up a practical method to consistently evolve three- and many-body interactions with the SRG 
transformation. In Ref. [1581, the first application of such a method to evolve three-body nuclear forces has been 
demonstrated, in calculations of A < 4 nuclei using a harmonic oscillator basis. 



3. Example: Two-dimensional Dirac-delta Contact Interaction 

In order to illustrate the general features of the SRG approach, in this section we discuss a simple example 
from quantum mechanics. We consider the SRG evolution of the hamiltonian for a system of two non-relativistic 
particles in two dimensions interacting via a Dirac-delta contact potential 1 159J - li,164J . A discretized version of the 
two-dimensional Dirac-delta contact potential was considered previously by Glazek and Wilson 116511 as a model to 
study the application of the SRG to an asymptotically free theory. Here we consider the continuous version. 

The Schrodinger equation in configuration space for the relative motion of two particles in two dimensions with 
an attractive Dirac-delta contact potential (in units such that h- 1) is given by 

- 'P(r) - Q-o 6^^\r) "Vir) = E ^(r) . (13) 

The coupling constant oq is dimensionless, since both the kinetic energy operator and the two-dimensional Dirac- 



delta contact potential scale as and so Eq. ( |13| l is scale invariant (i.e. there is no intrinsic energy or length 
scale). Although this equation can be easily solved in configuration space by regularizing the Dirac-delta contact 
potential with a distribution function, we choose to work in momentum space in order to estabhsh a more transparent 
connection with calculations implemented in the framework of SRG and EFT. 
The Schrodinger equation in momentum space is given by 



/ ^iP) f d^q m) = E (I)(p) , (14) 



where p is the relative momentum and 0(p) is the Fourier transform of the configuration space wave-function ^^(r). 



Scale invariance implies that if a bound-state solution exists (with £ < 0) then Eq. ( 14 1 will admit solutions for 
any negative energy, which correspond to a continuum of bound-states with energies extending down to -oo. The 
eigenvalue condition for the binding energy Eq - -E is given by 



27r Jo ip^ + Eo) 



The integral in the r.h.s. of Eq. ( 15 1 diverges logarithmically and so the problem is ill-defined. Nevertheless, using 



standard regularization and renormalization schemes we can obtain an exact analytic solution M159I . 



5.7. Standard Solution 



We start by regularizing the integral in Eq. {15\ using an ultraviolet momentum cutoff A, 

1 



27T Jo 



dp p 



(p^ + Eo) An 
Solving this equaion for the binding energy £0 we obtain 

£0 = {e'^"'"" - 1) 
6 



= — In 1 -I- — 



Eo 



(16) 



(17) 



Clearly, if the coupling constant q-q is fixed, then oo when A — > oo. To eliminate the divergence and produce 
a well-defined bound-state, we can renormalize the coupling constant by demanding that it depends on the cutoff A 
such that the binding energy £0 remains fixed as A — > 00, 

ao ^ ffA = — — T- . (18) 

Thus, the dimensionless renormalized running coupling q-a that characterizes the strength of the Dirac-delta contact 
interaction is replaced by a dimensionful parameter, the binding energy £0, which can be arbitrarily chosen and fixes 
the energy scale of the renormalized theory. This is a simple example of dimensional transmutation II166I . The original 
bare hamiltonian is scale invariant, but the renormalization process introduces an energy scale that characterizes the 
observables of the theory. It is also interesting to note that the renormalized theory is asymptotically free, since the 
renormalized running coupling vanishes in the limit A — > 00. 

We can use the renormalized interaction to evaluate other observables. The usual prescription for the calculations 
is to obtain the solutions with the cutoff in place and then take the limit A — > 00. If an exact calculation is implemented, 
the final results must be independent of the regularization and renormalization schemes. As an example, we evaluate 
the scattering phase-shifts from the solution of the LS for the T-matrix with the renormalized interaction, 

r(p,p';^2) = y(p,p')+ r d^q T(q,p';k^), (19) 

where k is the on-shell momentum in the center-of-mass frame. 

Only S-wave (I = 0) scattering occurs, since for the higher waves the centrifugal barrier completely screens the 
Dirac-delta contact potential. So, we can integrate over the angular variable to obtain 

1^ '{p,p -K-) = Vj^ '(_p,p)+ I dqq 

Jo 

where "'(/?,/?') = -a^/ln. The LS equation for the on-shell T-matrix (p - p' = k) is given by 

Tl=°\k, k;k')^-^-^ rf ">(fc, k; k') r dq q ^ \ . (21) 

' ' ' 2n In ^ ' Jo k^-q^ + ie 



Tl=\p,p'-e-) = Vt'\p,p') + dqq „^ Tl=%,p';k^) , (20) 

Jo k^ - q^ + I e 



Solving this equation, we obtain 



T'):'\k,k-k^) 



Inl^Uln^^^-^^" 



(22) 



K^ + k^ ^ 

The exact on-shell T-matrix is then obtained by taking the Umit A ^ 00: 

T^^='^\k, k; k^) = -2 [in [k^/Eo) - i nf^ . (23) 

A more convenient way to evaluate scattering observables is to use the /^-matrix (reactance matrix). The LS 
equation for the /^-matrix is similar to the one for the T-matrix, except that standing-wave boundary conditions are 
imposed for the Green's function. In this way, the ie prescription used in the LS equation for the T-matrix is replaced 
by the principal value, such that the /T-matrix is real. The S-wave LS equation for the /^-matrix is given by 



Kr(p,p';k')^Vr\p,p')^'P r dqq'^^^Kr%,p';k'), 

Jo 1^ 'J 



(24) 



where !P denotes the principal value. The relation between the exact on-shell A'-matrix and T-matrix is given by, 

K Hk,k;k^)^T-^ (k, k; k^) - i njl . (25) 
For simplicity, here and in what follows we drop the superscript 1 = 0. 
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Using the relation between the on-shell T-matrix and the phase-shifts, 



k cot 6(k) - ik 



2k 



T-\k,k-k^) 



(26) 



we obtain the exact phase-shifts. 



cot 6(k) 



-ln{k^/Eo) 



(27) 



It is important to note that keeping a finite momentum cutoff A in place leads to cutoff-dependent inverse power- 
law errors in the inverse on-shell T-matrix and hence in the phase-shifts. For A >> V^o, such errors can be 
evaluated perturbatively by double-expanding the difference between the A-dependent inverse on-shell T-matrix and 
the exact inverse on-shell T-matrix in powers of y/Eo/A and k/A, 



AT- 



T-'(k,k-k^)- 



T-\k,k-k^) 



(Eo + k^) (El-k') 



2A2 



4A4 



+ 



Ell_ 
A*-' A6 



(28) 



3.2. SRG Evolution 

The hamiltonian in momentum space for the relative motion of two non-relativistic particles in two dimensions 
with an attractive Dirac-delta contact potential can be written in the form 



^^(P,p') = riei(p,p') + v(p,p') 



(29) 



where TreKp, p') - /?^<5'^'(p - p') is the free hamiltonian (kinetic energy) and y(p, p') = -ai)/{27r)^ corresponds to 
the Fourier transform of the Dirac-delta contact potential in configuration space. One should note that the matrix 
elements of the Dirac-delta contact potential in momentum space are constants, clearly showing that the hamiltonian 
directly couples an infinite number of scales. Such couplings are precisely the source of the logarithmic divergences 
that appear in the Schrodinger equation, as shown in the previous subsection. 

We apply the SRG approach to renormalize this hamiltonian using Wegner's formulation. Choosing the generator 
= [Tre\, Hs] and integrating out the angular variable, we obtain the flow equation for the relative momentum space 
matrix elements of the potential in the basis of the free hamiltonian eigenstates, r,ei|p >= p^\p >, given by 



dVs(p,p') 
ds 



2s2 



--ip'-p") 



VAp, p) + I dqq(p^+ p'^ - 2q^)Vs(p, qWAq, p). 
Jo 



(30) 



For convenience, we consider here and in the rest of the paper that the flow parameter s is related to a similarity cutoff 
A with dimensions of momentum, such that s = A^'^. 

In principle, we could impose the boundary condition at sq = (/lo — » oo). 



HsMp^p') = H{p,p') = T,,Ap,p') + V(p,p') . 



(31) 



where T^^Ap, p') - p 5{p - p') and V(p, p') - -a(]l2n. However, the bare hamiltonian (i.e., with no similarity cutoff) 
produces logarithmic divergences and so the boundary condition should be imposed at some other value sq + 0, 
leading to dimensional transmutation. Formally, a renormalization prescription must be specified that allows one to 
obtain a well-defined hamiltonian at s - sq and thus fixes the underlying theory. 



A simple prescription is to use the renormalized running coupling aA derived in subsection (3.1 1 to define the 
potential in the bare hamiltonian. We introduce an ultraviolet momentum cutoff A and set 



y{p,p)^yK{p,p) 



' 2ti 



(32) 



Using this prescription, the bare hamiltonian becomes well-defined and the boundary condition can be imposed at 
So - 0. Note that the underlying theory is fixed at /io — > oo by fitting the binding energy Eq. 
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In previous works 0167II - II169I we solved the flow equation Eq. (30i perturbatively, using the idea of "coupling- 
coherence" M170I - M173I . We assumed a solution in the form of an expansion in powers of as,\l2n. 



n=2 



(33) 



where Ff'^{p,p') are irrelevant operators that vanish when p - p' - 0, and solved Eq. (30 1 iteratively, obtaining 
analytic expressions order-by-order for the perturbative approximations of the renormaUzed potential, such as: 



(a) marginal operator with initial coupUng a a. 



Vs,a(p,P )^-^ e 
2n 



(2) 

(b) marginal operator with second-order running coupling a\ j^, 

(2) 

2n 

(2) (2) 

(c) marginal operator plus second-order irrelevant operator F ' with second-order running coupling a i. 



(34) 



(35) 



(2) 



In 



In 



F^^A^P,P') 



-sip^-p'^f 



(36) 



Then, we used these approximations in a sequence of bound-state calculations and analyzed the scaling of the binding 
energy errors with the similarity cutoff' /I = 5 '^"^ at each order 

In Refs. 11671 1168 1 a comparison was made with bound-state calculations based on Lepage's treatment of EFT 
ifTSi . For simplicity, a separable effective potential was used that in momentum space consists of the Fourier transform 
of a Dirac-delta potential regularized by a smooth function of a momentum cutoff A and a series of approximately 
local effective interactions which correspond to the derivatives of the Dirac-delta potential. 



yfip,p') 



Co(A) + C2(A) 



2A2 



+ C4(A) 



(p'+p'l 
4A'* 



2 i2 
P P 



g-p-/2A2 g-/-/2A^ 



(37) 



Following the method described by Steele and Furnstahl 11741 11751. the parameters C,(A) were determined order- 
by-order by fitting the values for the exact inverse on-shell A'-matrix. The method consists in fitting the difference 
between the inverse on-shell K-maXiix obtained from the effective potential at a given order and the exact inverse 
on-shell A'-matrix to an interpolating polynomial in k^/A^, 



AK-^ = K,\k,k-k^)-K-\k,k-k^) = Ao 



A2 



■AA—r+- 

A4 



(38) 



and then minimizing the coefficients A, with respect to the variations of the parameters C,(A). 

The results from the error analysis have shown that for the EFT calculations the errors scale like inverse powers of 
the momentum cutoff A, as expected since the effective interactions are irrelevant operators. A systematic power-law 
improvement is obtained as each term is added to the effective potential and the corresponding parameter C,(A) is 
adjusted by the fitting procedure. For the SRG calculations, the error analysis displayed two different scaling regimes. 
When the similarity cutoff A is larger than the momentum cutoff A (introduced to define the bare hamiltonian) the 
scaling is similar to that for EFT, i.e. the errors are power-law. When the similarity cutoff A becomes smaller than 
the momentum cutoff A, in addition to the power-law errors there are inverse logarithmic errors that result from the 



truncation of the perturbative expansion Eq. ( 33 i for the renormaUzed potential at a given order in a, ^ 
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Although not trivial, the scaling behavior of the errors in the perturbative SRG can be qualitatively understood 
through fairly simple analytical calculations. By regrouping the terms in the SRG potential given by Eq. (33 i we can 
formally write it as a momentum expansion similar to the EFT potential given by Eq. ( [37] i, 



V.5,a(p,/?') 



(39) 



where, in this case, the expansion parameters gi(as /^) are analytic functions of the running coupling a, ^ 

gi(as,A) = fl;(«, A) as,A + bi(s. A) a^^ + ■ ■ ■ . 



(40) 



From this expansion, we can identify two interdependent sources of perturbative errors in the calculations of 
observables using the SRG potential truncated at a given order in the running coupling Us.a- First, there are inverse 
logarithmic errors introduced by the perturbative approximation for o-j a. Second, there are errors introduced by the 
perturbative expansion for the functions §,(Q'.s,a) powers of Os^a, which are a combination of inverse power-law 
and inverse logarithmic errors. Consider for instance the approximation (c) given by Eq. ([36l), with A and £0 fixed 

(1) (2.) 

(A >> Eq). The second-order running coupling a^^ and irrelevant operator F]^^ are given respectively by 

(2) _ Q'A 



l-^ [r + In (2iA4) - Ei (-2.?A4)] 



(41) 



and 



= ^[y + \n{2s p'^) - Ei(-2i p^ p'^)\ + ^ [r + ln(2iA'*) - Ei(-25A'*)] 

- \[y + \n{s [ip^ - K^f + ip'^ - A'f - ip' - p'^)']) 

-Ei [ip' - A'f + ip'^ - A'f - ip' - p'^f])] , 



(42) 



where Ei(x) is the exponential integral function. Expanding F*^' in powers of yfs p^ = p^/^^ and p'^ - p'^/A 



we obtain 



where 



vfiiP^p')- 



§0(«?i) + |2(«fi) 



2^2 



2/|4 



§o(«?i) 



(2) 

2n ' 

/„(2) ^2 i2 

87r2A2 



2a' 



,(2) 
,),A 



(43) 

(44) 
(45) 
(46) 



When A » A the second-order running coupling a^^\ runs little (approaching its asymptotic value qa) and so 



the functions g2{af\) and §4(0^^'^^). Evaluating the diff'erence between the truncated SRG potential V^j^^ip, p') and the 
initial potential Va(p,p') = -o:Al{2n) and expanding in powers of KjA, we obtain 



(2), 



w < ^-v(2), vr '^ (P' + P") 4aA (/,2p'2) 

SVaAp, p) = v\.(p,p)-Va(p,p)^--, 



8;r2^4 4^2^2 2^2 



2^4 



(47) 



One can clearly see that in this region perturbation theory in 6V,ia leads to errors that scale like powers of similar 
to EFT. In the limit /I — > 00 the results become exact. 

When A becomes smaller than A there is a crossover to a different scaling regime. The second-order running 
coupling a^^^j^ runs significantly in this region, and so the functions §2(0"^^^^) and §4(0^^^^, leading to errors that scale 



like a combination of powers of l/A and inverse powers of \n(A). 
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In this work we solve Eq.(30l numerically, obtaining an exact (non-perturbative) solution for the SRG evolved 
potential (apart from numerical errors). The relative momentum space is discretized on a grid of gaussian integration 
points (we have used 200 mesh points in our calculations), leading to a system of A^^ non-linear first-order coupled 
differential equations which is solved by using an adaptative fifth-order Runge-Kutta algorithm. 

In Figs. ([T]) and (j2]i we show respectively the contour and the surface plots for the evolution of the two-dimensional 
Dirac-delta contact potential with the similarity cutoff A = s^^^*, starting from an initial (bare) potential Va{p,p') 



defined by using the renormalized running coupling a a in Eq. ( 18 i fixed at a momentum cutoff A = 50 to give a 
binding energy £0 = 1 (in arbitrary units). The matrix elements of the initial potential Va=5{){p, p') are constants and 
we clearly see a systematic suppression of the off-diagonal matrix elements as the similarity cutoff A is reduced, such 
that the potential is driven towards a band-diagonal form. 
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Figure 1: (Color online) SRG evolution of the two-dimensional Dirac-delta potential for A = 50 and E(, = lin arbitrary units (contour plots). 




Figure 2: (Color online) SRG evolution of the two-dimensional Dirac-delta potential for A = 50 and £0 = 1 in arbitrary units (surface plots). 
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In Fig. |3]we show the results for the phase-shifts as a function of the laboratory energy ^lab calculated from the 
numerical solution of the LS equation for the /T-matrix with the initial two-dimensional Dirac-delta contact potential 
VA=5oip, p') and with the corresponding SRG potentials evolved to several values of the similarity cutoff A. 




Figure 3: (Color online) Phase-shifts for the initial two-dimensional Dirac-delta potential and the SRG potentials evolved to several values of A. 



In Fig. |4]we show the relative errors in the phase-shifts for the SRG evolved potentials with respect to the results 
for the initial potential Va=5o(p, p'), as a function of Zslab for several values of the similarity cutoff A (left) and as a 
function of the similarity cutoff A for several values of ^lab (right). Note that, in order to avoid mixing up with the 
errors from using a momentum cutoff A = 50 and with the errors from the numerical solution of the LS equation, 
we compare the phase-shifts for the SRG evolved potentials with those for the initial potential Va=5o(p, p'), rather 



than with the exact results given by Eq. ( 27 1. As one can observe, the phase-shifts obtained with the SRG evolved 
potentials are the same as those obtained with the initial potential, apart from relative differences smaller than 10 
due to residual numerical errors (which depend on the discretization and on the accuracy and relative error tolerance 
set in the Runge-Kutta solver used to evolve the potential). This result is expected, since the SRG transformation 
is unitary. One can also observe two different scaling regions in the error plots for fixed Zslab (right). When the 
similarity cutoff A is smaller than the momentum cutoff A(= 50) the errors are practically constant. When A becomes 
larger than A there is a crossover to a different scaling regime where the errors significantly decrease. 




10"' 10"' 10"' 10"' 10° io' 




Figure 4: (Color online) Relative errors in the phase-shifts for the SRG potentials evolved from the initial two-dimensional Dirac-delta potential. 
Left panel: as a function of £lab for several values of A; Right panel: as a function of A for several values of £lab- 

The change in the scaling behavior of the residual errors is a numerical artifact that can be understood in the 
following way. When A » A, the SRG potential evolves slowly and stays very near the initial potential. Only a 
few steps are required to integrate the differential equations using the Runge-Kutta solver, resulting in very small 
errors in the matrix-elements of the SRG evolved potential and hence in the phase-shifts. As A approaches A, the 
SRG potential starts to evolve more rapidly. An increasing number of steps is required to integrate the differential 
equations, resulting in larger errors in the phase-shifts. When A « A, the errors approach an asymptotic value. 
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4. SRG for the A^A^ Interaction 



In this section we review the main aspects involved in the application of the SRG approach to the A^A^ interaction 
II128141321 . We start from the hamiltonian in the center-of-mass frame for a system of two nucleons, which can be 
written in the form H - Trei + V, where Trei corresponds to the relative kinetic energy and V corresponds to the A^A^ 
potential. Here and in what follows we use units such that h - c - M - 1, where M is the nucleon mass. 

Using the generator z/, - [Trei, Hs\, Wegner's flow equation for the SRG evolution of the A^A^ potential is given by 



dVsip,p') 
ds 



-(p'-p'^fVAp,p' 



2 r 
) + - 

^ Jo 



dq ip^ + //2 - 2q^) Vsip, q) V^iq, p'). 



where Vs{p, p') is used as a brief notation for the projected A^A^ potential matrix elements, 

yf^^'^^^>(/.,/y) = ( p(LS)J;I\Vs\ p'{L'S)J-I) , 
in a partial-wave relative momentum space basis, | q(LS )J; I ), with normalization such that 

2 



1 = 



7!" Jo 



dqq^ \q(LS)J;I) {q(LS)J-I \. 



(48) 



(49) 



(50) 



The superscripts J, L{L'), S and / denote respectively the total angular momentum, the orbital angular momentum, the 
spin and the isospin quantum numbers of the A^A^ state. For non-coupled channels {L - L' - J), the matrix elements 
Vsip, p') are simply given by Vsip, p') = V\f''''^''\p, p'). For coupled channels (L, L' - J + 1), the Vsip, p') represent 
2x2 matrices of matrix elements for the diff'erent combinations of L and L': 



Vsip,p') = 



Vs ip, P ) 



Vs ip, P ) 



vr^^^^p,/) vr'''''\p,P')) 



(51) 



As a consequence of the choice of the SRG transformation generator rjs - {TrA, Hs], each interaction channel evolves 
with the cutoff A = '^^ independently of the other channels [I76|. 

The scattering observables for each A^A^ interaction channel can be calculated by iterating the corresponding SRG 
evolved potential through the LS equation for the partial-wave T-matrix (for simplicity, we drop the subscript s 
denoting the flow parameter and the superscripts denoting the quantum numbers of the A^A^ state): 



Tip,p';k^)^Vip,p')+- f dqq \/^'';'^\ Tiq,p';k^). 

^ Jo k-^ - q^ + 1 £ 

The phase-shifts are then obtained from the relation Sik^) = 1 - 2ik Tik,k; k^). For the non-coupled channels, 

Sik^)^t^^i2i6'f), 



and for the coupled channels (using the Stapp parametrization II177I ). 

cos (2e') exp(2 / 5'j_^) i sin(2e') exp(/ 6'j_^ + i 

i sin (2e') exp(/ + i cos (2e') exp(2 / 



Sik^) 



(52) 



(53) 



(54) 



where the superscripts s and t stand respectively for the spin singlet states and spin triplet states, J - L + S is the total 
angular momentum of the A^A^ state and e' is the mixing parameter. 

For convenience, in our calculations we iterate the potential through the LS equation for the partial-wave A'-matrix 



Kip,p';k^)^Vip,p')+-'P r dqq^^^^Kiq,p'-k^), 

^ Jo k^- q 



JO - 

The relation between the T-matrix and the /T-matrix on-shell is given by: 

R-^k, k; k^) = T-\k, k;k^)-ik. 
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(55) 



(56) 



In order to illustrate the application of the SRG to the A^A^ interaction, we solve Eq. (48 i numerically to evolve the 
Nijmegen potential |55 1 in the '5o channel. In Figs, (jsj and ^ we show respectively the contour and the surface plots 
for the SRG evolution of the Nijmegen potential in the 'S'o channel. As one can observe, the initial potential (A — > oo) 
exhibits strong matrix elements outside the diagonal, extending up to high-momenta (large energy differences). As 
the potential evolves with the similarity cutoff A, the off-diagonal matrix elements are systematically suppressed while 
the high-momentum components are gradually concentrated near the diagonal and the low-momentum components 
(attractive) are enhanced. Therefore, the Nijmegen potential evolved through the SRG tends to a band-diagonal form. 





Figure 6: (Color online) SRG evolution of the Nijmegen potential in the So channel (surface plots). 
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In Fig. [Tjwe show the A^A^ phase-shifts in the Sq channel as a function of the laboratory energy ^lab calculated 
from the numerical solution of the LS equation for the partial-wave K-maXiix with the initial Nijmegen potential and 
with the corresponding SRG potentials evolved to several values of the similarity cutoff A. As expected, the phase- 
shifts obtained with the SRG evolved potentials are the same as those obtained with the initial potential, apart from 
relative differences smaller than 10"' due to numerical errors. This result is obtained for all energies, even those not 
constrained to scattering data (larger than /slab ~ 35QMeV). 
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Figure 7: (Color online) Phase-shifts in the So channel for the initial Nijmegen potential and the SRG potentials evolved to several values of A. 



Now, we test the decoupling of low-energy observables from high-energy degrees of freedom for the Nijmegen 
potential evolved through the SRG, following the analysis introduced by Bogner et al. 1 1291 fT30J . Such an analysis 
consists in applying an exponential regularizing function to the potential, which suppresses the contributions to the 
observables from matrix elements Vsip, p') with p, p' larger than a given momentum cut A;tnax, 



Vf--'"Hp,p') = exp[-(pVCxn Vsip,p') exp[-(p'VCx)"] ^ 



(57) 



where n is an integer number that specifies the smoothness of the regularizing function. 

We apply the exponential regularizing function to cut the initial Nijmegen potential (A — > oo) in the ^So channel 
and the corresponding SRG potential evolved to a similarity cutoff A - 2.0 fm ' and then calculate the phase-shifts 
as a function of the laboratory energy /slab- In Fig- |8]we show the results obtained by cutting the potentials at 



^max - 2.2 fm ' (left) and k^-ix = 1.0 fm ' (right), with n 
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Figure 8: (Color online) Phase-shifts in the So channel for the initial Nijmegen potential and the corresponding SRG potential evolved to a 
similai'ity cutoff A = 2.0 fm"' with the application of an exponential regularizing function. Left panel: using n= 8 and /t,,,^.,. = 2.2 fm"'; Right 
panel: using n=8 and k^ax = 1.0 fm"' . 
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As one can observe, the phase-shifts obtained for the cut initial potential are completely modified when compared 
to those obtained for the uncut potential. For the SRG evolved potential cut at kranx - 2.2 fm ' (above A), the 
phase-shifts at low energies agree with those obtained for the uncut SRG evolved potential and deviations start to 
occur near ^lab ~ 2 k^^^/M. When the SRG evolved potential is cut at k^^^ -1.0 fm"' (below A), there is no 
agreement between the phase-shifts practically at all energies. One should note that the phase-shifts obtained for the 
cut potentials (both the initial and the SRG evolved) tend to zero for energies /slab > 2 k^^^/M. Therefore, an explicit 
decoupling between the low- and high-momentum components is verified for the SRG evolved potential, such that the 
high-momentum matrix elements above the similarity cutoff A can be removed without causing a significant distortion 
in the calculations of low-energy observables. This result is similar to those described in Refs. ||129|[r30J for other 
NN potentials. It is important to emphasize that different decoupling patterns are expected by using the generator 
T/j = [Gs, Hs] with other choices for the operator Gj, such as Gj - diag(Hs). 



5. The Subtracted Kernel Method 



5.7. Formalism 

We begin the description of the SKM approach II 134111 391 by considering the formal LS equation for the J-matrix 
of a two-body system, which can be written in operator form as 



T{E) = V + V G^{E) TiE) = V [l + G,;(£) T{E)\ 
= y + T(E) g;(E) V^[l+ T(E) Gt,(E)] V , 



(58) 



where V is the interaction potential and Gq{E) is the free Green's function for the two-body system with outgoing- 
wave boundary conditions, given in terms of the free hamiltonian Hq by 



G+(£) = (£ - Hq + ie)-' 



(59) 



For singular potentials, like those containing Dirac-delta contact interactions and their derivatives, Eq. ( [58] l becomes 
ill-defined due to the ultraviolet divergencies that appear in the implicit momentum integration. In the SKM approach 
a regularized and renormalized scattering equation is derived by performing subtractions in propagator at a certain 
energy scale. 

We first consider a singular potential containing a regular term and a Dirac-delta contact interaction, in which case 
only one subtraction is enough to get a finite scattering amplitude. In momentum space, this potential can be written 
in the form V = Vreg + Co, where Co is the strength of the contact interaction. Using Eq. (58 1, the potential V can be 
formally written in terms of the T-matrix at a given energy scale -/u^ (in units such that h - c - 2m - 1, where m is 
the reduced mass of the two-body system): 



V = r( V) [i + G,;(V) T{-t?)Y = [i + T{-^?) G+( V)]"' n-^?) . 



(60) 



For convenience we choose a negative energy for the subtraction scale, such that the free Green's function G^(-fr) is 
real (this is not a restriction, since the method works as well for a positive energy subtraction scale). 

Replacing the potential V in Eq. (58 1 by its expression in terms of T{-fi^) given in Eq. (60 1, we obtain 



T(E) = r( V) [l + Go(-M^) r( V)]"' + T(E) Gt,(E) T{-i?) [l + Gj( V) r( V)]" 



(61) 



Now, multiplying both sides of Eq. (61 



the subtracted kernel LS equation forme T-matrix, 



from the right by |^1 -i- Gq (-//^) r(-/z^)j and re-arranging the terms, we obtain 



T{E) = T{-^?) + T{-^?) [GliE) - Gl{-^^)\ T{E) , 



(62) 



which has the sam e op erator structure as the Eq. (58 1 for the formal LS equation, but with the potential V replaced 
according to Eq. ( 60 1 by the T-matrix at the subtraction scale -/i^ and the original free propagator replaced by a 
propagator with one subtraction at the same scale. 
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The subtracted kernel LS equation provides a finite solution for the T-matrix at any given energy E, once its 
value at the subtraction scale -yu^ is known. Thus, the input for the solution of Eq. (62 1 is T(-fp-), which is called 
"driving term" and contains the physical information apparently lost due to the removal of the propagation through 
intermediate states at the scale -fi^. A simple ansatz consists in considering that the driving term is given by 

T(-f/) = V(-f/) = Keg + Co(-ju') , (63) 

where Co{-fi^) is the renormalized strength of the contact interaction at the subtraction scale -fj.^, which is fixed by 
fitting experimental data for scattering observables and, therefore, encodes the physical information. 

In order to clarify the application of the SKM scheme through a simple example, we consider the scattering of two 
nucleons with the LO "pionless" EFT potential, which consists of just a pure Dirac-delta contact interaction, V - Co, 
and gives an analytic solution for the T-matrix. Using the ansatz T{-fi^) = V(-//^) = Co(-//^), the subtracted kernel 
LS equation for the T-matrix is given in operator form by 

T(E) = Co(-M^) + Co(-f/) [g;(E) - G,;(V)] T(E) . (64) 

As pointed out in subsection due to the centrifugal barrier only S-wave scattering can occur for a pure Dirac- 
delta contact interaction. Using a partial-wave relative momentum space basis with normalization given by Eq. (jSOjl, 
we obtain the subtracted kernel LS equation for the matrix elements of the T-matrix in the S-wave channel 



k^ - q^ + ie 



- q- 



r''(q,p';k^). 



(65) 



where k - Ve is the on-shell momentum in the center-of-mass frame and the superscripts s and t stand respectively 
for the spin-singlet channel ^Sq and the spin-triplet channel i. 



Solving Eq. ( [65] l by matrix inversion, we obtain 

^•'(p,p';k^)^ 



(fi + i k) 



(66) 



One should note that in the case of a pure Dirac-delta contact interaction the T-matrix is momentum independent. 

The renormalized strength of the contact interaction, C^'i-ji'-), can then be fixed by using the scattering length 
flj , as a physical input, which is defined by the relation a, , = r'''(0, 0, k^ - 0). From Eq. (66 1 we obtain 



1 



(67) 



which leads to a yu independent J-matrix, given by 

r-'ip,p;k^)^ 



— + ik 



(68) 



Such a result is equivalent to the one obtained from LO "pionless" EFT calculations using a momentum cutoff A to 
regularize the divergent integral in the formal LS equation and taking the limit A — > oo [^3"64"6T '69', '701 . Other 
regularization schemes, such as dimensional regularization with power divergence subtraction (PDS) |56, 56, 58 1 and 
momentum subtractions Il68ll69ll74ll75l . also lead to the same result in LO. Thus, this example clearly demonstrates 
that, in the case of a potential containing just a pure Dirac-delta contact interaction, a finite T-matrix can be obtained 
by solving the subtracted kernel LS equation with only one subtraction. 

For a general singular potential containing Dirac-delta contact interactions plus regular momentum dependent 
terms and/or derivative contact interactions, like the A^A^ potential to a given order in ChEFT, the subtracted kernel LS 
equation must be solved numerically. In this case, a finite T-matrix can still be obtained, but using the fitting procedure 
described before to fix the renormalized strengths in the driving term Ti-fP') generates a residual dependence with 
respect to p in the T-matrix, and hence in the scattering observables. As we will show in subsection (5.2i, a more 



rigorous procedure can be used in order to obtain scattering observables that are strictly invariant under the change of 
the subtraction point, which is based on a renormalization group equatio n for the running of the driving term T{-p^) 
with the subtraction scale -yU^. Furthermore, as we will see in subsection (5.3 i, when the original potential specifically 
contains derivative contact interactions the SKM scheme must be generalized to allow for multiple subtractions, 
leading to a highly non-trivial dependence of the driving term on the subtraction scale 
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5.2. Renormalized Potential and Fixed-point Hamiltonian 

Once the renormalized strength of the contact interaction C()(-yU^) in Eq. ( 63 1 is fixed, and so the driving term 
Ti-n^) is known, a renormalized potential Vn can be formally defined from Eq. (60 1 fl361: 



(69) 



Multiplying both sides of Eq. (69 1 from the left by [l + T(-fi^) Gq(-iu^)^ and re-arranging the terms, we obtain 
the integral equation (in operator form) 



which formally relates the renormalized potential V?j to the driving term T{-fi^). 



(70) 



Replacing V by V-r in Eq. ( 58 i, we obtain the LS equation for the renormalized T-matrix: 



Tk{E) = Vk + Vk G^(£) Tk . (71) 
The renormalized potential Vn is not well-defined for singular interactions. Nevertheless, for a driving term T{-fi^) 



like the one defined by Eq. ( 63 1, containing a regular term plus a Dirac-delta contact interaction, Eq. ( 7 1 1 gives a finite 
solution for the renormalized T-matrix, T^iE), which is equivalent to the one obtained from the subtracted kernel LS 



equation Eq.(62i, i.e. Tn{E) - T(E). This result provides a po^sfenon justification for the formal manipulations of 
Vk and V used to obtain Eqs. (|62]) and ( fTO] ). In numerical calculations it is more convenient to regularize the integral 
implicit in Eq. ( 70 1 by introducing an ultraviolet momentum cutoff A, such that we can work with a well-defined 



renormalized potential VniA). Using this potential, we construct the LS equation for a cutoff regularized T-matrix, 



Tk(E; A) = Vk(A) + y«(A) (5+(£; A) Tk(E; A) , 



(72) 



where Gq(E; A) denote the Green's function in the regularized integral. 

In the limit A — > oo the result obtained by solving Eq. (72 1 with the cutoff regularized potential VKiA.) should be 
the same as the one obtained from the direct solution of Eq. (71 1 with Vn- The same procedure could be applied for 
the numerical solution of the subtracted kernel LS equation Eq. (62 1. It is important to emphasize that the momentum 
cutoff A used in this context just plays the role of an instrumental regulator for the numerical integration. The relevant 
scale parameter is the subtraction scale at which the physical information is introduced through the renormalized 
strength of the contact interaction Coi-fi^) in the driving term T^-ji^). 

The subtraction scale -//^ is arbitrary, and so all scattering observables should not depend on its choice. In order 
to fulfill this condition, the T-matrix must be independent of the subtraction scale -yu^: 



dT 



= . 



From Eqs. (71 1 and (73 1 we find that the renormalized potential V-r is also independent of -fi 



= 0, 



(73) 



(74) 



which implies that the renormalized hamiltonian Hn - Hq + V'r corresponds to a renormalization group fixed-point, 
i.e. it is stationary with respect to the change of the subtraction scale -fi^. 



Using Eqs. (70i and (74i or Eqs. (62 1 and (73 i, a renormalization group equation can be obtained for the driving 
term T{-^p-) 1 136 1 in the form of a non-relativistic Callan-Symanzik equation (NRCS) Ill78til80l . 



do 



()(-^') 2, 
^^-^ ^ ' 



(75) 



with the boundary condition, imposed at some reference scale /i, given by )|^^^ = Vi-jl ) - Vreg + C()(-jj ). 

The NRCS equation provides a definite rule for the evolution of the driving term T{-fi^) with the sliding scale yU, from 
the reference scale //, such that the T-matrix remains invariant under the dislocation of the subtraction point. 
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5.3. Recursiveness of the SKM for Multiple Subtractions 

Using a notation that will be convenient when considering multiple subtractions, the LS equation for the T-matrix 
with one subtraction, Eq. ( 62 1, can be written as 

T{E) = y^'>( V) + V) -I?) T(E) , (76) 

where y'^'(-/z^) = T{—n^) is the driving term and G|(£; —fJ^) = Gq{E) - G^i-jj.^) denote the subtracted Green's 
function with one subtraction, which can also be written in the form 

G|(£; -m') = (-yu' - E) Gti-f?) Gl{E) = F,{E- -f?) G+(£) , (77) 

The function Fi (E; —fj^) can be regarded as a form factor that modifies the free Green's function Gq{E) by introducing 
a factor proportional to q^^ in the momentum integration. Such a form factor is enough to regularize the divergence 
generated in the LS equation by a singular potential containing just a pure Dirac-delta contact interaction. However, for 
singular potentials containing terms that generate higher-order power divergences in the LS equation, like derivative 
contact interactions, the SKM scheme must be generalized. As shown in Refs. 1135111371 . in this case it is necessary 
to perform multiple subtractions in the kernel of the LS equation by using an iterative procedure. 

For a general number of subtractions n, we define a n-fold subtracted kernel LS equation, given by 

T{E) = V^"\E; V) + V^"\E; -//) G^(E; -//) T{E) . (78) 

The n-fold subtracted Green's function G^(E; -j-P') is defined by 

G:(E; -ti^) = [(-//2 - E) G^(-//2)]" G^E) = F„(E; -f?) G^E) . (79) 

Note that the «-fold form factor F„(E; -ju^) introduces a factor proportional to q^^" in the momentum integration, 
thus regularizing divergences to order q^"^K The driving term V^'"\E; -fi^) for the «-fold subtracted LS equation is 
recursively constructed through an iterative procedure, starting from y*''(-//^). The recursion formula is given by 

y('"'(£; -t?) = y("'>(£; -t?) + vl^^^i-/?) , (80) 

where 

V^'"\E; -^?) = [l - (-n^ - ET-' V^"'-^\E- -fi^) G^i-i^Yl V^"'-'\E- , (81) 

and the term y''"^^(-^J^) contains the higher-order singular interactions that generate divergent integrals which can 
be regularized by performing m subtractions. One should note from Eq. (80i the that the driving term y''"'(£'; -jJ^) 
at each iteration is derived in two steps. First, we calculate V^'"\E\ -yU^) from V'^'"^^\E; -jJ^), solving the integral 



equation obtained by manipulating Eq. (81 



y(""(£; -I?) = y("'-i'(£; -f?) + (-^i^ - £)'""' V"'-^\E; -fi^) Gli-i?)"' V^'"\E\ -i?) . (82) 

Then, we introduce the corresponding higher-order singular interactions in the driving term by adding ys™g(-A'^). 

As in the case of one subtraction, once the driving term V^"\E; -ji^) with n subtractions is known we can formally 
manipulate Eq. (78 i to define a renormalized potential, 

Vk = [\ + V^"\E; -fi^) (g,;(£) - G,t(£; -fi^))]'' y(">(£; -fi^) , (83) 



which is iterated in the LS equation Eq. (71 1 for TuiE) and leads to a renormalization group fixed-point hamiltonian. 
In the same way, a NRCS equation for the driving term y*"'(£'; -jj.^) can be obtained from the invariance of the n-fold 
subtracted J-matrix with respect to the subtraction scale -/i^, 

with the boundary condition given by V^"\E; -fJ.^)\fi^p = V^"\E; -p^) imposed at some reference scale Jl. 
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As an example to illustrate how the SKM scheme works when multiple subtractions are performed, we consider the 
scattering of two nucleons in the S -wave channel with the NLO "pionless" EFT potential, which consists of a Dirac- 
delta contact interaction plus a second-order derivative contact interaction. In a partial-wave relative momentum space 
basis, the matrix elements of such a potential are given by: 



(85) 



Even thought the NLO potential contains only terms of order q^, when iterated in the LS equation it generates integrals 
that diverge as much as q^. Therefore, we have to construct a 3-fold subtracted kernel LS equation in order to get a 
finite result for the T-matrix. The iterative procedure starts from the T-matrix with one subtraction at -/i^, such that: 



V''\p,p';-fi')^Co(-i^') 



The next step is to derive the driving term V^^K We calculate V'^' from V*^'' through Eq. ( 82 1: 



2 r°° 

y<2)(p,p';fc2;_^2)^y(Vp';_^2)+_ dq q 

^ Jo 

Solving this equation by matrix inversion, we obtain 

V^^\p,p';k^;-f/)^ 

with 



^2 + ^2x v^^Hp,q;-^i^) 



y<2)(^,p';fc2;_^2) 



1 



+ /o 



^ Jo 



dq q 



(fj.^ + ^^2)2 2jU 

Since there are no higher-order singular terms to be added at this iteration, we simply have 

r 1 /..2 , ;,2^^-I 



1 ^(p^ + k^) 



Co(-P^) 



2p 



(86) 



(87) 



(88) 



(89) 



(90) 



The last step of the iterative procedure is to derive the driving term V^^\ First, we calculate V*-** from V^^^ through 
the integral equation 



2 r°° 

V^^'>(p, p'; k^; -p^) = y(2>(p, p'- k^- -p^) + - dqq^ 

^ Jo 

which is solved by matrix inversion giving 



p^+k^\ V''^Hp,q;k^;-p^) 



p^ + ^2 



with 



V^^\p,p';k^;-p^) = Co(k^;-p^)^ 

2 r°° 

h=h{k^-p^)^- dqq 
^ Jo 



1 



Coi-p^) 



(-A*' - q^) 



+ k + h 



V^\q,p'-k\-p^), (91) 



(92) 



2 {p^ + k^f _{p^ + k^f 



Ou2 + qiy (2pf 



Then, we add the term vf„g(p,/7'; -p'^) = Cii-p^) (p^ + p"') to obtain 



y(3>(p,p';^2._^2) 



1 



_^ (//2 + k^) _^ip^+ k^Y 



C2(V) (p'+p ) 



[Coi-p^) 2p {2pf 
The resulting 3-fold subtracted LS equation for the T-matrix is then given by: 

(p^ +k^^ V'^^Hp,q;k^;-p^) 



T(p,p';k') = V^'\p,p';k'--p') + 



2 r 

^ Jo 



dq q 
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p^ + q^ 



k^ - q^ + ie 



T{q,p';k^) 



(93) 



(94) 



(95) 



Using the method described in Refs. 



analytically. The driving term V'^* given in Eq. ( 94 1 can be written in the form of a two-term separable potential, 

I 

i,j=0 



, suitable for separable potentials, Eq. ([95]l can be solved 



(96) 



where Aij(k , // ) are the matrix elements of 



Co C2 
Ci 



with Co = Co(k^; -fi^) and C2 = Cii-fi^). 

The solution of Eq. ( 95 1 can then be written in the form 



(97) 



T{p,p'-k^) = P^' ^7(^';/"') ■ 

/,>=o 



The unknown matrix Tik , ^ ) satisfies the equation 



T{k^;i?) = Hk';n') + A(^^;/) J{k'■^i') T(k';fi') , 



r,2. ,,2- 



2. .,2\ ,,2x rrn.2. ,,2\ 



where 



with 



J„ = Uk^;^?) 



dq 



iP- + k^ \^ ^2«+2 



(« = 0,1,2). 



(98) 

(99) 
(100) 

(101) 



\//^ + / k^ - + ie 

Note that the factor F^ik^; -jj^) = (//^ + k^)^ /(jj^ + q^)^, introduced by the subtraction procedure, regularizes the 
otherwise-divergent integrals J„{k^;fp-) for n = 0, 1 , 2. 



Solving Eq. (|99|l for T{k^; fi^) and substituting in Eq. (|98| yields 

T(p,p';k^) 
Using the identity 



Co + C2 (p^ + p'^) + Cl [72 - (/ + //') Ji + -/q] 
\ - Cq Jq -2 C2 J\ - C^Ui Jq - 



2 r 

^ Jo 



dq f(q) 



k^-q^ 



2n+2 2 r°° fl^""*"^ 

— = -;p ^?/(?)t| — J -ik^"^\ 

+ le n Jq - g'^ 



and the relation (valid for n > 0) 



2 

^ Jo 



dq f(q) 



q 



2n+2 



k^-q^ 

the integrals 7„ can be written in the form 



2 r 

^ Jo 



dq f(q) q^" + k^ - P dq f(q) ^' 
^ Jo 



2n 



k^-q^ 



Ja^WQ + k^W-ik\ 
j^^Wx+k^WQ+k'^W- i P ■ 
J2^W2 + k^ Wi+k'^Wo + k^W- i k^ , 



with 



W„^W,Ak\p')^-- rdql^^^] q'"; W^W{1^■,^?)^-'P r dq(^ 
n Jo \p^ + q^j n Jq + 



2 , ^2x3 



q^j k^-q^ 



(102) 



(103) 



(104) 



(105) 
(106) 
(107) 



(108) 
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The inverse on-shell 3-fold subtracted T-matrix (p - p' — k) can then be written in the form 



T{k, k- k^) Co + C2 W2 + k^ C2 [2 -C2W1] 



- (Wo + k^W) + ik, 



(109) 



One should note that the unitarity of the 3-fold subtracted r-matrix is ensured by the term ; k. 

The renormalized strengths Co(-fi^) and Cii-ji'^) can be fixed by fitting the experimental values of the scattering 
length a and the effective range r^. Expanding Eq. \\Q9) in powers of k yields 



1 



= Ao[Co, C2; -ti'] + Ai [Co, C2; V] k' + (3(r) + ik. 



T{k,k-M) 

with the coefficients Ao[Co, C2; -^P'] and A2[Co, C2, -^^^] given by 



Ao[Co,C2;-ju1 



2 64 -H 64 Co ^l+\6C2^i^ + 10 Co C2 /^^ - 8 C^ - 5 Co C5 //^ 



64 Co - 24 Cl A<5 - 15 Co 



(110) 



(111) 



Ai[Co,C2; V] 



256 [32 C2 - 8 C^ + d 1/ + 2C\ ju'] 

[64 Co - 24 Cl - 15 Co jjfif 
64 Co [112 C2 yu - 82 C2 / H- 5 C^ + 10 C^ yu'O] 

[64 Co - 24 Cl - 15 Co //S]^ 
4 C2 [l28 C2 ^i - 632 C2 + 25 C^ + 50 C\ 
[64 Co - 24 C2 - 15 Co a/«]^ 



Matching Eq. ( 1 10 1 to the first two terms in the effective range expansion (ERE), 



1 



• [k cot 5{k) - i k] 



-- + - r,k^ + 0(k'^)-ik 
a 2 



we obtain a system of two coupled non-linear equations for Cq{-p^) and C2(-A(^), 



Ao[Co,C2; V] = 



^ Ai[Co,C2; V] = ■ 



(112) 



(113) 



(114) 



(115) 



The solutions of these equations for the renormalized strengths Co(-fP') and C2(-/U^) are highly non-trivial functions 
of the scale yu. In Fig. Q we show one of the solutions obtained by using parameters for the A^A^ scattering in the 
'5o channel as physical input, i.e. an unnaturally large scattering length a - -23.7 fm and a positive effective range 
rg - 2.7 fm. For comparison, we also show the result obtained for Co(-/i^) in the calculation with just a pure Dirac- 
delta contact interaction, given by Eq. (67i. As one can observe, for values of the scale /i up to //max ~ 0.8 fm"' 
the solutions for Co(-/U^) and C2(-yU^) are real. For /i > /i,nax the solutions become complex, corresponding to a 
non-Hermitian driving term, and thus lead to the violation of unitarity. This result implies that in order to fit a positive 
effective range while preserving the unitarity of the subtracted T-matrix we cannot use an arbitrarily large value for the 
scale fi. Such a behavior is similar to what is found in NLO "pionless" EFT calculations using the cutoff regularization 
scheme ||59l l6T|j64l l67l 16911701 and can be related to the Wigner bound on the effective range [1811, a general result 
based on the physical principles of causality and unitarity which states that for a Hermitian potential of range R the 
effective range is constrained by 



r,<2 



R' 

R -(- — 



(116) 



If we take R ~ l/fi, then for a - -23.7 fm and 
with the results of our calculations. 



2.7 fm we obtain the condition jj. < 0.78 fm which is consistent 
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|i (MeV) H (MeV) 

Figure 9: (Color online) Renormalized strengths Coi-fJ^) and C2(-/i^) adjusted to fit the '5o channel scattering length and effective range. 




k (MeV) k / 1^ 

Figure 10: Left panel: (Color online) Left panel: Phase-shifts in the '5o channel as a function of the on-shell momentum for the LO and NLO 
potentials renormalized through the SKM procedure, compared to the results from the Nijmegen PWA and from the ERE to order k-; Right panel: 
relative errors in the phase-shifts (with respect to the results from the ERE) as a function of (k/fi) for the LO and NLO potentials (/i = 130 MeV). 



In the left panel of Fig. ( 10 1 we show the A^A^ phase-shifts in the '5o channel obtained for the LO potential from 
Eq. ( [68| ) and those obtained for the NLO potential from Eq. ( |109| l with several values of the scale fi, compared to the 
results from the Nijmegen partial-wave analysis (PWA). We also show the A^A^ phase-shifts obtained from the ERE 
to order k^, which provides a very good fit to the results from the Nijmegen PWA up to on-shell momenta k of order 
m„ ~ 140 MeV. As one can observe, the phase-shifts for the NLO potential renormalized through the SKM procedure 
are in good agreement with those obtained from the ERE up to a given value of the on-shell momentum that becomes 
larger as the scale yU increases. Note that for very low on-shell momenta (up to ~ 40 MeV) the phase-shifts are nearly 
independent of the scale /i. Such results would be expected, since the renormalized strengths Co(-/i^) and C2(-/U^) in 
the NLO potential are fixed by matching the ERE to order k- and so dependence on the scale fi in the inverse on-shell 
3-fold subtracted T-matrix should start at order (k/jj)'^. This can be clearly verified from the log-log plots for the 
absolute value of the relative errors in the phase-shifts with respect to the results from the ERE, shown in the right 
panel of Fig. ( 10 1, whose slopes are given by the dominant power of k/fi in the errors. 

It is important to emphasize that once a reference scale fi is specified and the renormalized strengths are fixed to fit 
the observables used as physical input (e.g., the scattering length and the effective range), the subtraction point -fi^ can 
be changed without modifying the results for the calculated phase-shifts, as long as the driving term V^^\p, p';k^; -jJ^) 
in the 3-fold subtracted LS equation for the J-matrix is evolved through the NRCS equation Eq. ( [84| with the 
boundary condition given by V^^Hp, p'; k^\ -fi^)\fi^p = V^^\p, p'; k^; -fP'). 
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We have just shown two simple analytical examples of how to renormalize A^A^ singular interactions using the 
SKM approach with one and multiple subtractions. We close this section with some remarks about the advantages 
and disadvantages of the SKM formalism and its possible extension to renormalize 3A^ interactions. 

As pointed out in the introduction, an advantage of the SKM formalism is that it can be extended to any derivative 
order of the contact interactions, through an interative process involving multiple subtractions. Such an iterative 
process allows for a systematic treatment of ChEFT A^A^ potentials up to higher-orders, which include interactions 
that generate higher divergences. This also makes the SKM convenient to test alternative power-counting schemes 
in which higher-order contact interactions are promoted to lower-order [801. Furthermore, the renormalization group 
invariance of the SKM, ensured by the evolution of the driving term through a NRCS-type equation, allows one to 
arbitrarily move the renormalization scale (i.e. the subtraction point //) while maintaining the observables invariant. 
One can even slide the renormalization scale beyond the limits imposed by the Wigner bound, which is not possible 
within the cutoff regularization scheme. A disadvantage of the SKM is that the computational load increases as 
more subtractions are performed, since the recursive process requires more matrix inversions. So, while the cutoff 
regularization scheme demands only one inversion to compute the T-matrix, the SKM with n subtractions requires 
n - \ inversions to compute the driving term V*"' plus one inversion to compute the «-fold subtracted T-matrix. 

The SKM formalism presented here provides a powerful method to renormalize two-body interactions. Similar 
subtractive renormalization methods 11143141461 . developed in the context of EFT for three-body systems with contact 
interactions, suggest that the SKM formalism can be generalized to treat three-body interactions. As in the SKM 
approach with one subtraction, in the methods described in Refs. M143H1461 renormalized amphtudes are obtained by 
performing a subtraction in the three-body scattering equations, so we believe that an extension of the SKM approach 
to derive renormalized ChEFT potentials including 3A^ interactions might work as well. Once obtained, such SKM 
renormalized 3A^ potentials can then be used as the starting point for the SRG evolution. Investigating the possibihty 
of applying the SKM formalism to the 3A^ system is certainly a project we intend to pursue in forthcoming works. 



6. SRG Evolution of A^A^ Potentials in the SKM Approach 



6.1. SKM for the NN Interaction in LO ChEFT 

In this subsection we present a detailed and systematic analysis of the SKM approach applied to the A^A^ interaction 
in LO ChEFT, which consists of the OPEP plus a Dirac-delta contact interaction. This is a simple and convenient 
example to perform a detailed investigation, since only one subtraction is required to renormalize the LO potential. 
Such an application has already been described in Ref. [134.1 for the singlet '5o channel and the triplet \ Di 
coupled channel. Here, we restrict our calculations to the '5o channel only. In Ref. II 1371 the SKM approach was 
applied to the A^A^ interaction in the '^o channel with the inclusion of a second-order derivative contact interaction in 
addition to the LO ChEFT interaction, requiring the use of a 3-fold subtracted kernel LS equation in order to regularize 
the divergences. The implementation of the method with the full A^A^ interaction in NLO ChEFT, which includes 
the NLO two-pion exchange potential (TPEP), was considered in Ref. [11381 in the case of the '^o channel. In Ref. 
|fT39 [ the method was applied to describe peripheral waves in uncoupled channels with the inclusion of the NNLO 
TPEP, requiring the use of a 4-fold subtracted kernel LS equation. 

The A^A^ interaction in LO ChEFT is given by 



V(p,p') 



d;-(p' -^&2-{p' -P) , 1 



(p' -fPi^ + ml 



1 - r 1 ■ T2 



+ C' 



3 + Tl ■ T2 



(117) 



where and tr, (i - 1 , 2) are the usual isospin and spin Pauli operators, ga - \ .25 is the axial coupling constant, 
f„ - 93 MeV is the pion weak-decay constant and m^j - 138 MeV is the pion mass. The coefficients Cg and Cy 
correspond to the strengths of the contact interactions respectively for the '5o and \ Di channels. 



In a partial-wave relative momentum space basis with normalization given by Eq. ( 50 1, the matrix elements of the 
LO ChEFT potential projected in the '5o channel are given by 



V(p,p')^Vi„{p,p') + Cl 



(118) 
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where Vi„(p,p') are the corresponding OPEP matrix elements, given by 

VUp,p') = ^(2- fdx- ^^^^^ -]. (119) 

^^^fn \ J -I p'- + p'^- 2pp'x + ml I 

We implement the SKM procedure by using the /T-matrix instead of the T-matrix. We also introduce an ultraviolet 
momentum cutoff A which is convenient for the numerical calculations when we further consider the evolution of the 
SKM renormalized potential through the SRG transformation. Such a cutoff allows one to work with a finite-range 
gaussian grid of momentum integration points, which is required for an efficient numerical solution of Wegner's flow 
equation using a Runge-Kutta solver. As pointed before, the cutoff A must be regarded as an instrumental regulator 
for the numerical integrations, whose effects on the calculated quantities should vanish in the limit A — > 00. 

The subtracted kernel LS equation for the cutoff regularized /T-matrix in the '^o channel, with the momentum 
cutoff A included through a step function 6{A - q), is given by 

Kip,p';k^;A)^mp,p'--p^-A)+-r f dqq^e{K-q) f4^) '^''\^'~i'^^ K(q,p'-,k\K) . (120) 

7T Jo +q I k^-q^ 

Note that in the limit A — > 00 the result obtained by solving this equation should be the same as the one obtained by 
solving the equation for the unregularized /T-matrix, i.e. K(p,p';k^;A-^ 00) ^ K(p,p';k^). 



The driving term V^^\p, p'; -ju^; A) at the subtraction scale -fj.^ in Eq. ( 120 1 is defined by the ansatz 

y('>(/,,p'; V; A) = K(p,p';-f/;A) = ViJp,p') + A) ■ (121) 

The renormalized strength of the contact interaction Cj)(-jU^; A) is fixe d at t he scale yu for a given momentum cutoff 
A by fitting the scattering length a, in the '^o channel. We solve Eq. ( 120 1 numerically and evaluate the scattering 



length from the on-shell /T-matrix (p - p' - k) at E - k - through the relation as - K(Q, 0; k - 0; A), adjusting 
the value of Q(-yU^; A) so as to reproduce the experimental value - -Ti.l fm with a given accuracy. 

In Figs. ((ll|i and ((l2ll we show the scattering length a, as a function of the strength Cy(-yL/^; A), respectively for 
several values of the scale with the momentum cutoff A fixed and for several values of the momentum cutoff A 
with the scale ji fixed. As one can observe from both plots, the position of the resonance pole (associated with the 
Feshbach resonance 018211 due to presence of a shallow virtual bound-state in the '^o channel) converges to a limit 
value as the scale ^ and the momentum cutoff A increase. The same happens with the renormalized strength of the 
contact interaction (the value adjusted to reproduce Oj = -23.7 fm), which is very close to the position of the pole. 



In Fig. (lob we show the running of the renormalized strength Cq{-ij^\ A) with the scale yu for several values of the 
momentum cutoff A (left) and with the momentum cutoff A for several values of the scale ji (right). For a given value 
of the momentum cutoff A the dependence of (^^(-yU^; A) on the scale yU vanishes asymptotically in the limit yu ^ 00 
and so the SKM scheme leads to the same results as would be obtained by using the (sharp) cutoff regularization 
scheme. Similarly, for a given value of the scale yu the dependence of Cq{-^'^\ A) on the momentum cutoff A vanishes 
asymptotically in the limit A — > 00, i.e. Cy(-yU^; A ^ 00) — > Cy(-yU^), with Cj)(-yu^) absorbing the divergences and 
leading to the correct value of the scattering length a ,. This result clearly shows that the role played by the momentum 
cutoff A is purely instrumental, since in the limit A —> 00 all ultraviolet divergences in the LS equation are regularized 
by the subtraction procedure itself, providing a finite /T-matrix. 
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Figure 1 1 : (Color online) Scattering length Uj as a function of for several values of the scale /j with the momentum cutoff A fixed. 
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Figure 12: (Color online) Scattering length cij as a function of for several values of the momentum cutoff A with the scale p fixed. 
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Figure 13: (Color online) Renormalized strength of the contact interaction C^^-j?", K). Left panel: as a function of the scale /i for several values of 
the momentum cutoff A; Right panel: as a function of the momentum cutoff A for several values of the scale /i. 
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Although simple, the numerical procedure used in the SKM scheme to fit the scattering length Os requires a very 
delicate tuning of the strength C^i-fi^; A). As shown in Figs. (Ill and ( 12 1, the resonance width becomes smaller 



as the scale /i and the cutoff A increase (up to the point where the position of the pole stabilizes), such that a higher 
tuning precision is necessary to achieve a given accuracy. 

A) must be tuned, defined as 



In Fig. 



14 



we show the results obtained for the relative precision to which C()(-jU' 
AC()/Cq, in order to fit the scattering length in the '5o channel with an accuracy of 1%, as a function of the scale jj 
for several values of the cutoff A. We also show the results obtained for each cutoff A in the limit /i ^ oo, which 
correspond to using the (sharp) cutoff regularization scheme. 

As expected, for a given cutoff A the required precision becomes higher as the scale fi increases (corresponding to 
a smaller value of ACq/Cq), asymptotically approaching the precision required when using the cutoff regularization 
scheme. One should also note that in the limit A — > oo (no cutoff) the tuning of the strength Cp(-jt/^; A) should become 
numerically impossible for scales // larger then some maximum value. 
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Figure 14: (Color online) Relative precision ACg/C^ to which the strength Cg(-/j-; A) must be tuned in order to fit the scattering length in the 'So 
channel with an accuracy of 1%. The solid lines with symbols correspond to the SKM results as a function of the scale ^ for several values of the 
cutoft' A. The horizontal dashed lines correspond to the values for each cutoft' A obtained by using the cutoft' regularization scheme (;j — > oo). 



We now turn to the calculation of the renormalized potential Vii{p,p'; A). Once the renormalized strength of the 
contact interaction C^^i-fi^; A) is fixed at the scale ju for a given momentum cutoff A, we can obtain the renormalized 
potential Vk{p, p'; A) from the driving term V^^\p, p'; -p^; A) by solving an integral equation which is derived from 



Eq. (70 1 (replacing the T-matrix by the /T-matrix), 



Vn{p,p';A)^Y 



,7(1) 



ip,p';-p^-h)--'P 

7T 



dq q 6(A - q) 



W\p,q--p^;A) - 



Vn{q,p'\A) , 



(122) 



Note that we keep the instrumental cutoff A in place, so as to obtain a SKM renormalized potential which will be 
numerically more convenient to evolve through the SRG transformation. 

In principle, the renormalized potential VKip, p' ; A) obtained from the solution of Eq. 



122 1 should be independent 



of the subtraction scale -jj , since by construction it corresponds to a renormalization group fixed-point operator. 
However, a residual dependence on the scale p is generated by the fitting procedure used to fix the renormalized 



strength of the contact interaction C^i-p , A). As pointed out in subsection (5.2 1, in order to obtain a renormalized 
potential which is rigorously invariant with respect to the change of the subtraction scale -p^ the driving term must be 



evolved via the NRCS equation, with the boundary condition V^^\p, p'\ -p^; A)|^_ 
some reference scale p where the renormalized strength C'^{-p^; A) is fixed. 



y*'^(/7, p'; -p^; A) specified at 
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In Fig. (15 1 we show the driving terms V^^\p,Q; -fp-; A) (left) and the corresponding renormalized potentials 



-1 



Vk(p, 0; -jj. , A) (right) at p' — 0, for several values of the scale yu with the momentum cutoff fixed at A = 25 fm 
At each value of p., the driving term was obtained by tuning the renormalized strength C'^i-p^; A) to fit the scattering 
length fl^ in the ^Sq channel with the same accuracy, corresponding to an absolute error of 10"^. In Fig. (16 1 we show 
the relative differences between the driving term at a scale p and at /i — > oo (left), 

AV"^ ^ [WKp, 0; V; A) - WKp, 0; -ex,; A)] 
yd) y(i)(;',0;-oo;A) 

and between the renormalized potential at a scale p and at — > oo (right), 

AVk ^ [V^ip, 0; -ph A) - Vk(p, 0; -oo; A)] 
Vk ~ yR(p, 0; -oo; A) 



(124) 



as a function of the momentum p for several values of the scale p (top panels) and as a function of the scale p for 
several values of the momentum p (bottom panels). As expected, the driving term V'-^Xp, 0; -p^; A) changes with the 
scale p in such a way that the renormalized potential V^ip, 0; -p^; A) remains approximately invariant. In the limit 
— > oo, the driving term matches the renormalized potential, both becoming independent of the scale p. One can 
also observe that for a fixed momentum p the relative differences scale approximately as l/p^ in the range of scales p 
considered in the calculations, both for the driving term and the renormalized potential. 

Once the driving term V^^\p, p'; -p^; A) is fixed at the subtraction scale -p^ for a given momentum cutoff A, 
we can obtain the A^A^ phase-shifts for the LO ChEFT potential from the numerical solution of the subtracted kernel 



LS equation for the partial-wave /T-matrix, Eq. ( 120 1. In Fig. ( 17i we show the results for the A^A^ phase-shifts 
in the '^o channel as a function of the laboratory energy Elab calculated for several values of the scale p with the 
momentum cutoff A fixed at 25 fm"' (left) and for several values of the momentum cutoff A with the scale p fixed 
at 25 fm"' (right). The results obtained for the phase-shifts from the subtracted kernel LS equation strongly deviates 
from the Nijmegen PWA, similar to the well known results obtained for the LO ChEFT interaction in the '^o channel 
renormalized through other schemes. 



In Fig. ( 18 I we show the relative differences between the phase-shifts in the channel calculated at a scale p 



and at yu — » oo (with the momentum cutoff A fixed at 25 fm '), 

AS[p] ^ [^(£lab; -a/'; A) - ^(£lab; A)] 
S[p] 5(£'lab; -oo; A) 



(125) 



as a function of /slab for several values of p (left) and as a function of p for several values of /slab (right)- In the range 
of scales p considered in the calculations, the relative differences A6[p]/d[p] scale approximately as £'lab(°^ k^) for 
fixed p and as l/p^ for fixed /slab- In the limit yu — > oo, the phase-shifts should become independent of p, matching 
the results obtained by using the cutoff regularization scheme. 



In Fig. il9\ we show the relative differences between the phase-shifts in the '5o channel calculated at a momentum 
cutoff A and at A — > oo (with the scale p fixed at 25 fm"'), 

AS[A] ^ [^(£lab; -a^'; A) - ^(£lab; ^)] ^^26) 

S[A] (J(£'lab;-/U^;oo) 

as a function of /^lab for several values of A (left) and as a function of A for several values of /^lab (right). In the 
range of momentum cutoffs A considered in the calculations, the relative differences A5[A]/(5[A] scale approximately 
as £'lab('^ k^) for fixed A, similar to the scaling of A6[p]/d[p]. On the other hand, for fixed ^lab they scale as 1/A. 
As expected, in the limit A — > oo the phase-shifts should become independent of A. It is important to remind that 
the cutoff A is just an instrumental regulator, whose effects on the phase-shifts could be eliminated by using standard 
methods to take the limit A — > oo, such as the Richardson extrapolation [1831. 

We emphasize again that by evolving the driving term V^^\p,p'; -p^;A) from a reference scale p to the scale p 
via the NRCS equation, instead of using the fitting procedure to fix the renormalized strength C^^(-p^; A), the residual 
dependence of the renormalized potential V'Rip, p'; ; -p^; A) and the phase-shifts on the scale p should vanish. In this 
case, a non-trivial dependence of the driving term on the scale p is expected in order to ensure the invariance of the 
subtracted /T-matrix, which would lead to a more compUcated scaUng behavior. 
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Figure 15: (Color online) Driving term V'-'Mp, 0; -/i^;A) (left) and renormalized potential Vk(p, 0; -/J"; A) (right) for several values of the scale fi 
with the momentum cutoff fixed at A = 25 fm"' . 
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Figure 16: (Color online) Relative differences AV^'V^*'' for the driving term (left) and AV«/Vk for the renormalized potential (right). Top panels: 
as a function of the momentum p for several values of /j; Bottom panels: as a function of p for several values of the momentum p. 
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Figure 17: (Color online) Phase-shifts in the '5o channel as a function of the laboratory energy Slab for the LO ChEFT potential evaluated from 
the subtracted kernel LS equation for the A'-matrix. Left panel: for several values of the scale /i with the momentum cutoff A fixed at 25 fm"'; 
Right panel: for several values of the momentum cutoff A with the scale fi fixed at 25 fm"' . 




Figure 18: (Color online) Relative dilferences AS[ji]/S[ji] (with the momentum cutoff A fixed at 25 fm ). Left panel: as a function of Slab for 
several values of the scale /i; Right panel: as a function of the scale for several values of Elab ■ 




Figure 19: (Color online) Relative differences Ai5[A]/5[A] (with the scale fi fixed at 25 fm ). Left panel: as a function of Z?lab for several values 
of the momentum cutoff A; Right panel: as a function of the momentum cutoff A for several values of £lab • 
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6.2. SRG Evolution of the SKM-LO ChEFT potential 

When considering the application of the SRG approach to evolve the effective A^A^ potential in LO ChEFT, we use 
as an initial potential the fixed-point renormalized potential Vk derived by implementing the SKM scheme for the LO 
ChEFT interaction. We solve Eq. (48 1 numerically, with the boundary condition set at 5 = (/I — > oo) such that the 
initial potential Vs=o{p,p') is given by Vii{p,p') (as pointed before, for numerical convenience the calculations are 
performed by using an ultraviolet momentum cutoff A that acts just as an instrumental regulator). 
In Figs. 



(20 1 and (21 



we show respectively the contour and the surface plots for the SRG evolution of the 
SKM-LO ChEFT potential in the '5o channel. The initial potential Vii(p,p') was calculated by fixing Coi-jj^) at 
p. - 25 fm ' and using a momentum cutoff A = 25 fm '. As one can observe, the initial SKM-LO ChEFT potential 
has non-zero off-diagonal matrix elements which extend up to high momenta. As the similarity cutoff A is lowered, 
the off-diagonal matrix elements are systematically suppressed while the low-momentum components are enhanced, 
similar to what happens for the Nijmegen potential. Thus, the result clearly shows that the SKM-LO ChEFT potential 
evolved through the SRG transformation is driven towards a band-diagonal form. 
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Figure 20: (Color online) SRG evolution of the SKM-LO ChEFT potential in the So channel (contour plots). 




In Fig. (22 1 we show the A^A^ phase-shifts in the '^o channel as a function of ^lab calculated from the numerical 
solution of the (formal) LS equation for the partial-wave /iT-matrix with the initial SKM-LO ChEFT potential and 
the corresponding potentials evolved through the SRG transformation to several values of the similarity cutoff A. As 
expected, the SRG evolved potentials yield the same results as the initial potential for all energies, apart from relative 
differences smaller than 10"' due to numerical errors. 
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Figure 22: (Color online) Phase-shifts in the 'So channel as a function of the laboratory energy Elab for the initial SKM-LO ChEFT potential and 
the corresponding SRG potentials evolved to several values of A. 



In Fig. ( 23 I we show the relative errors in the phase-shifts for the SRG evolved potentials (with respect to the 
results for the initial potential) as a function of Zslab for several values of the similarity cutoff A (left) and as a 
function of the similarity cutoff A for several values of ^lab (right). Similar to what is obtained for the SRG evolution 
of the two-dimensional Dirac-delta potential, two different scaling regions can be observed in the error plots for fixed 
Elab (right), with a crossover at /I ~ A = 25 fm"^ 





Figure 23: (Color online) Relative errors in the phase-shifts for the SRG potentials evolved from the initial SKM-LO ChEFT potential. Left panel: 
as a function of £lab for several values of A; Right panel: as a function of A for several values of i^LAB- 



In Fig. ( [24| we show a comparison between the momentum-space matrix-elements of the Nijmegen and the 
SKM-LO ChEFT potentials in the '^o channel evolved through the SRG transformation. As one can observe, the 
potentials remain clearly distinct when the similarity cutoff is lowered to /} = 1 fm Such a behavior is expected, 
since the two initial potentials are not phase-shift equivalent (as shown before, the A^A^ phase-shifts obtained for the 
SKM-LO ChEFT potential strongly deviates from those obtained for the Nijmegen potential). Nevertheless, one can 
see a tendency of the matrix-elements of the two SRG evolved potentials to become similar at low-momentum, in 
particular the off-diagonal matrix-elements. 
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Figure 24: (Color online) SRG evolution of the Nijmegen and the SKM-LO ChEFT potentials in the 'So channel. Left: Diagonal matrix-elements; 
Right: Off-diagonal matrix-elements. 



Since the description of low-energy A^A^ observables provided by ChEFT potentials becomes systematically more 
accurate (i.e., approach the Nijmegen PWA results) at higher-orders in the chiral expansion (e.g., NLO, NNLO and 
N-^LO), low-momentum universality as described by Bogner et. al. 0128| [T33I is expected to emerge as higher-order 
contributions are included in the initial SKM ChEFT potential, i.e. the low-momentum parts of the Nijmegen and the 
SKM ChEFT potentials evolved through the SRG transformation should flow to quantitatively similar forms as the 
similarity cutoff is lowered. This must be verified by explicit calculations of SRG evolved higher-order SKM ChEFT 
potentials, which we will perform in a future work. 

In order to analyze the decoupling between low- and high-momentum components in the SRG evolved SKM-LO 
ChEFT potential, we use the exponential regularizing function defined by Eq. (57 1. In Fig. (25 i we show the results 
for the A^A^ phase-shifts in the '5o channel as a function of /slab obtained by cutting the initial SKM-LO ChEFT 
potential and the corresponding SRG potential evolved to a similarity cutoff X-2 fm"' at k-^a^ = 3.5 fm"' (left) and 
^max = 1-8 fm"' (right), with n - %. For the initial SKM-LO ChEFT potential, the phase-shifts obtained for the cut 
potential are completely modified in comparison to those obtained for the uncut potential. 

The qualitative decoupling pattern for the SRG potential evolved from the SKM-LO ChEFT potential is similar 
to that obtained for the SRG potential evolved from the Nijmegen potential. The phase-shifts remain practically 
unchanged at low energies when the SRG evolved potential is cut above X and significantly change for all energies 
when it is cut below X. 
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Figure 25: (Color online) Phase-shifts in the 'So channel for the initial SKM-LO ChEFT potential and the con'esponding SRG potential evolved 
to /I = 2 fm"' with the application of an exponential regularizing function. Left: n=8 and = 3.5 fm"'; Right: n=8 and = L8 fm"'. 



A more quantitative analysis of the decoupling pattern can be performed by evaluating the relative errors in the 
phase-shifts for the cut SRG evolved potential (with respect to the results for the corresponding uncut potential) as a 
function of the momentum cut fcmax, which provide a measurement of the degree of decoupling 1 130|. In Fig. pS] ) we 
show the log-log plots for the relative errors obtained for the cut SRG potentials evolved from the SKM-LO ChEFT 
potential (left panels) and from the Nijmegen potential (right panels), respectively for several values of ^lab. the 
regularizing function parameter n and the similarity cutoff A. As one can observe, there are three different regions 
in all error plots, both for the SRG potentials evolved from the SKM-LO ChEFT and from the Nijmegen potential. 
As expected, when femax < VA/£'lab/2 the relative errors go to one, since in this case the matrix-elements of the 
cut potentials vanish such that the corresponding phase-shifts go to zero. When k^^,, > A the relative errors scale 
as a power-law (l/^max)^". indicating a perturbative decoupling regime. This power-law scaling behavior is clearly 
evidenced by the straight lines with slope -2n in the error plots. It is important to note that the crossover to the 
power-law scaling region typically occurs at fc^ax slightly above A, independently of the values of /slab and n. In the 
intermediate region the scaling behavior is not definite, becoming different for the SRG potentials evolved from the 
SKM-LO ChEFT and from the Nijmegen potential. 

In Fig. (27 1 we show the relative errors in the phase-shifts for the SRG evolved potentials cut at k,^^,^ - 10 fm"' 
with n = 8 as a function of ^lab for fixed A - 2 fm"' (left) and as a function of A for fixed ^lab = 100 MeV (right). 
As one can observe in the left panel, the dependence of the relative errors on ^lab is weak for Z^lab < 100 MeV. In the 
right panel one can observe that the dominant relative errors scale as a power-law A^" in the range 2 fm"' < A < 6 fm"' 
for the SRG potentials evolved from the SKM-LO ChEFT and in the range 2 fm"' < A < 5 fm"' for the SRG potentials 
evolved from the Nijmegen potential. 

The general results we have obtained in the quantitative analysis of the decoupling pattern for the SRG potentials 
evolved from the SKM-LO ChEFT potential and from the Nijmegen potential in the '^o channel are similar to those 
described in Ref 1 130| for the SRG potentials evolved from the chiral N^^LO potential with a 500 MeV cutoff ll42l . 
However, there is a difference which is worth to mention. As shown in Ref. [113011 . the change in the phase-shift 
relative errors for the cut SRG evolved potentials saturates when the similarity cutoff A approaches the underlying 
500 MeV (~ 2.54 fm"') momentum cutoff of the initial chiral N^LO potential, in particular the change in the position 
of the shoulder in the error plots indicating the crossover to the power-law scaling region. Here, as one can observe in 
the bottom panels of Fig. (26 1, there is much less saturation for similarity cutoffs /i up to ~ 5 fm"', because the initial 
SKM-LO ChEFT potential has no underlying cutoff (except for the large instrumental cutoff A = 25 fm"'), going to a 
small constant value for momenta larger than ~ 6 fm"', and the initial Nijmegen potential has an effective momentum 
cutoff ~ 5 fm"' (which results from the combination of the several gaussian cutoff parameters used to regularize the 
terms associated to the different meson exchange processes included in the potential 1551 ). Indeed, as one can see in 
the right panel of Fig. ( 27 1, the saturation effect is observed only for values of A larger than ~ 5 - 6 fm"' . 
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Figure 26: (Color online) Relative errors in the phase-shifts for the cut SRG evolved potentials as a function of the momentum cut, ^max, for several 
values of Slab, 'i and A. Left: SKM-LO ChEFT potential; Right: Nijmegen potential. 
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Figure 27: (Color online) Relative errors in the phase-shifts for the SRG evolved potentials cut at /:max = 10 fm with n = 8. Left: as a function 
of -Elab for fixed A = 2 fm"' ; Right: as a function of A for fixed -Elab = 100 MeV. 



7. Summary and Conclusions 

So far, we have performed a detailed investigation of the similarity renormalization group (SRG) evolution of the 
A^A^ interaction in leading order chiral effective field theory (ChEFT) within the subtractive kernel method (SKM) 
renormalization approach. In the following, we summarize the results of our study and present our main conclusions. 

First, we considered the two-dimensional Dirac-delta contact interaction, aiming to illustrate the application of 
the SRG through a simple example and to fine-tune the parameters of the Runge-Kutta solver used in our code to 
obtain an accurate numerical solution of Wegner's SRG flow equation. The matrix elements of the two-dimensional 
Dirac-delta potential in momentum space are constants, thus generating ultraviolet logarithmic divergences when 
iterated in the LS or the Schrodinger equation. Nevertheless, exact analytic solutions can be obtained by introducing a 
momentum cutoff and renormalizing the coupling constant that characterizes the strength of the contact interaction in 
order to remove cutoff dependence in the observables, leading to dimensional transmutation and asymptotic freedom. 
We solved Wegner's flow equation for the two-dimensional Dirac-delta potential re normalized by demanding that the 
binding-energy remains fixed as the momentum cutoff is removed and verified that the SRG evolved potential is driven 
to a band-diagonal form as the similarity cutoff A is reduced. From the numerical solution of the LS equation for the 
^T-matrix we calculated the phase-shifts as a function of Zslab for the initial two-dimensional Dirac-delta potential 
and for the corresponding SRG potentials evolved to several values of the similarity cutoff A and verified the unitarity 
of the SRG transformation up to relative numerical errors smaller than 10"'. Two scaling regions were observed in 
the error plots for fixed ^lab: the relative errors are nearly constant for A < A and significantly decrease for A > A. 
Such a change in the scaling behavior of the errors is a numerical artifact related to the number of steps required by 
the Runge-Kutta solver to evolve the potential. 

Next, we considered the SRG evolution of the high-precision Nijmegen A^A^ potential in the ^Sq channel. We also 
verified the band-diagonalization of the SRG evolved potential and, by evaluating the phase-shifts, the unitarity of 
the SRG transformation (up to relative numerical errors smaller than 10"'). Then, following the method described by 
Bogner et al. M129II130 T we used an exponential regularizing function to smoothly cut the SRG evolved potentials 
above a given momentum k^^^. We verified that, as a consequence of the decoupling of low-energy observables from 
high-energy degrees of freedom, the phase-shifts calculated with the SRG evolved potentials cut at a momentum 
above the similarity cutoff A agree with those obtained for the corresponding uncut potentials up to /slab ~ 2 k^^^/M. 

Then, we presented a discussion of the SKM renormalization approach. As a first example, we considered the 
A^A^ interaction in LO pionless EFT, which simply consists of a pure Dirac-delta contact term. In this case, only one 
subtraction in enough to renormalize the interaction and the subtracted kernel LS equation can be solved analytically 
to obtain the T-matrix. The renormalized strength of the contact interaction is fixed by fitting the experimental value 
of the scattering length. Our result for the scattering amplitude in the '^o channel, obtained from the subtracted kernel 
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LS equation, is independent of the subtraction scale and exactly matches the one obtained by solving the formal 
LS equation regularized using a momentum cutoff A and then taking the limit when such a cutoff goes to infinity. 
Next, as an example to illustrate the application of the SKM renormalization approach when multiple subtractions 
are performed in the kernel of the LS equation, we considered the A^A^ interaction in NLO pionless EFT. In this case 
three subtractions are required, since the NLO interaction also includes a first-derivative contact term, and so a 3-fold 
subtracted kernel LS equation must be recursively constructed. Using the method described in Refs. 1641 |69l l87l . 
the 3-fold subtracted kernel LS equation can still be solved analytically to obtain the scattering amplitude. The 
renormalized strengths of the two contact interactions are determined by matching the calculated scattering amplitude 
and the ERE to order (where k is the on-shell momentum), such that their values at a given subtraction scale fi are 
fixed by fitting the scattering length and the eff'ective range in '5o channel. We found that in order to get real values 
for the renormalized strengths, thus preserving the unitarity of the T-matrix, the value of subtraction scale fi must 
be less than ~ 0.8 fm a result which is consistent with the Wigner bound on the eff'ective range. Moreover, the 
dominant relative errors in the phase-shifts for the NLO potential with respect to the results obtained from the ERE to 
order k^ scale, as expected, like {kljSf. 

We then addressed the SRG evolution of chiral A^A^ potentials renormalized within the framework of the SKM. In 
our calculations, we implemented the SKM procedure by using the /T-matrix instead of the T-matrix. For numerical 
convenience when evolving the potential through the SRG transformation, we introduced an instrumental momentum 
cutoff" A in order to work with a finite-range gaussian grid of integration points. 

First, we performed a detailed analysis of the renormalization of the A^A^ interaction in LO ChEFT through the 
SKM approach. The LO ChEFT potential consists of the OPEP plus a Dirac -delta contact interaction and so only one 
subtraction is required. From the solution of the subtracted kernel LS equation in the '^o channel, we evaluated the 
on-shell /T-matrix 2LiE - k^ -Q and by fitting the scattering length in the '^o channel we investigated the running of 
the renormalized strength of the contact interaction with the subtraction scale ju and with the instrumental cutoff" A. 
For a fixed A the renormalized strength increases with ^ and in the limit — > oo it converges asymptotically to the 
same value that is obtained by using the sharp cutoff regularization scheme. Similarly, for a fixed yu the renormalized 
strength increases with A and in the limit A — > oo its dependence on A vanishes asymptotically, evidencing the 
instrumental nature of such a regulator. We also found that the precision to which the renormalized strength must be 
tuned in order to fit the scattering length with a given accuracy becomes higher as both ji and A increase. 

Next, we evaluated the driving term and the corresponding fixed-point renormalized potential for several values 
of the subtraction scale ^ with the instrumental cutoff A fixed. As expected, the driving term runs with ^ such that 
the fixed-point renormalized potential remains approximately invariant. In the limit yu — > oo the driving term becomes 
independent of //, matching the fixed-point renormalized potential. We found that the relative differences between 
the driving term at a given n and at ju — » oo scale like for fixed momenta. The relative diff'erences for the 
corresponding fixed-point renormalized potentials are about two orders of magnitude smaller, also scaling like Ij^P'. 
The residual dependence of the fixed-point renormalized potential on // is a consequence of the fitting procedure used 
to fix the renormalized strength of the contact interaction and should be eliminated by evolving the driving term from 
a reference subtraction scale to the subtraction scale // via the NRCS equation. Then, we calculated the phase-shifts in 
the '5o channel as a function of /slab for several values of the subtraction scale /i and the instrumental cutoff' A. The 
relative differences between the phase-shifts calculated at a given ju and at — > oo (with A fixed) scale like £'lab('^ k^) 
for fixed ju and like for fixed iiLAB- The relative diff'erences between the phase-shifts calculated at a given A and 
at A ^ oo (with jd fixed) scale like £'lab(<^ k^) for fixed A and like 1/A for fixed /Jlab- 

Finally, we moved to the calculations of the SRG evolution of the SKM-LO ChEFT potential in the '^o channel. 
We solved Wegner's flow equation using as an input the fixed-point renormalized potential derived through the SKM 
approach. As expected, the SRG evolution systematically suppresses the off-diagonal matrix elements as the similarity 
cutoff" A is lowered, driving the potential towards a band diagonal form. By evaluating the phase-shifts, we also 
verified the unitarity of the SRG transformation (apart from relative numerical errors smaller than 10"^). Like in the 
case of the SRG evolution of the two-dimensional Dirac-delta potential, two scaling regions were observed in the 
relative error plots for fixed /slab, with a crossover at /I ~ A. By comparing the matrix-elements of the Nijmegen 
and the SKM-LO ChEFT potentials evolved through the SRG transformation, we verified that the potentials remain 
distinct as the similarity cutoff A is lowered, showing only a tendency to become similar at low-momentum. This 
result is expected, since the two initial potentials are not phase-shift equivalent. Nevertheless, universality is expected 
to emerge as higher-order contributions in the chiral expansion are included in the initial SKM ChEFT potential. 
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As previously described for the Nijmegen potential, we analyzed the decoupling of low-energy observables from 
high-energy degrees of freedom by using an exponential regularizing function to cut the potential above a given 
momentum ^max- We found a qualitative decoupling pattern for the SRG evolved SKM-LO ChEFT potential similar 
to that obtained for the SRG evolved Nijmegen potential. In order to perform a quantitative analysis, we evaluated 
the relative errors in the phase-shifts for the cut SRG evolved potentials as a function of fcmax- From the log-log error 
plots we verified that for fc^ax > ^ the relative errors scale like a power-law (1/fcmax)^" both for the SKM-LO ChEFT 
and the Nijmegen potentials, indicating a perturbative decoupling regime. 

The main purpose of this work was to set up a basis for a comparative study of ChEFT A^A^ potentials using the SRG 
transformation, particularly through the SRG decoupling pattern analysis introduced by Bogner et al. ||129|[T30i . The 
simple examples considered here indicate that the SRG provides a powerful tool for analyzing the scale dependence 
of effective A^A^ interactions and for estimating the uncertainties in calculations of A^A^ systems. 

Our next steps, to be pursued in forthcoming works, will be to study the SRG evolution of ChEFT A^A^ potentials 
renormalized via the SKM approach up to higher-orders in the chiral expansion (NLO, NNLO and N^LO) and in other 
partial-wave channels (including calculations of phase-shifts and deuteron observables) and to consider the effects of 
using different generators for the SRG transformation. We will compare the results for the SKM ChEFT A^A^ potentials 
with those for ChEFT A^A^ potentials renormalized using cutoff regularization and investigate the interplay between 
the power counting (and hence the systematic improvement) and the non-perturbative renormalization using these two 
schemes. Another important issue to be addressed in such a study will be to verify how low-momentum universality 
of the SRG evolved potentials emerges as the higher-order contributions in the chiral expansion are included in the 
initial potentials and what are the effects of using different regularization schemes. We will also investigate the 
renormalization group invariance in the SKM formalism, by solving the non-relativistic Callan-Symanzik equation 
(NRCS) for the driving term of the subtracted LS equations (both in pionless EFT and ChEFT), and its possible 
extension to renormalize 3A^ interactions. Then, we intend to perform a detailed comparative analysis of several 
accurate ChEFT A^A^ potentials, including the chiral N^LO potentials of Epelbaum, Glockle and MeiBner I.39J and of 
Entem and Machleidt Il42l . the chiral N-'LO potential renormalized using the SKM approach and the brazilian chiral 
potential |53|, as well as the Argonne V18 potential (54] and the Nijmegen potential fSF] which are to be used as a 
baseline for the analysis. One of our motivations in this study is to investigate the effects on A^A^ observables due to 
the mid-range correlated two-pion exchange terms present in the brazilian chiral potential, which we expect may be 
better understood through the scale-dependence analysis provided by the SRG. 
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